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Abstract In broad astrophysical contexts of large-scale 
CN gravitational collapses and outflows and as a basis for 
various further astrophysical applications, we formulate 
^ and investigate a theoretical problem of self-similar mag- 
netohydrodynamics (MHD) for a non-rotating polytropic 
gas of quasi-spherical symmetry permeated by a com- 
C\| pletely random magnetic field. Within this framework, 
we derive two coupled nonlinear MHD ordinary differen- 
' ' tial equations (ODEs), examine properties of the mag- 
^ ^ nctosonic critical curve, obtain various asymptotic and 
I global semi-complete similarity MHD solutions, and qual- 
2 ify the applicability of our results. Unique to a magne- 
tized gas cloud, a novel asymptotic MHD solution for 
^ a collapsing core is established. Physically, the similar- 
I— I ity MHD inflow towards the central dense core proceeds 
in characteristic manners before the gas material even- 
' ^' tually encounters a strong radiating MHD shock upon 
Q,^ impact onto the central compact object. Sufficiently far 
1^ away from the central core region enshrouded by such an 
MHD shock, we derive regular asymptotic behaviours. 
CO We study asymptotic solution behaviours in the vicinity 
\^ of the magnctosonic; c;ritic;al cnirvc and cictcrminc smooth 
MHD eigensolutions across this curve. Numerically, we 
f-^ — construct global semi-complete similarity MHD solutions 
that cross the magnctosonic critical curve zero, one, and 
^ two times. For comparison, counterpart solutions in the 
• ^ case of an isothermal unmagnetized and magnetized gas 
flows are demonstrated in the present MHD framework 
H at nearly isothermal and weakly magnetized conditions. 



For a polytropic index 7 = 1.25 or a strong magnetic 
field, different solution behaviours emerge. With a strong 
magnetic field, there exist semi-complete similarity solu- 
tions crossing the magnctosonic critical curve only once, 
and the MHD counterpart of expansion-wave collapse so- 
lution disappears. Also in the polytropic case of 7 = 1.25, 
we no longer observe the trend in the speed-density phase 
diagram of finding infinitely many matches to establish 
global MHD solutions that cross the magnctosonic crit- 
ical curve twice. 

Keywords magnetohydrodynamics • planetary neb- 
ulae: general • stars: AGB and post-AGB • stars: 
formation • stars: winds, outflows • supernovac: general 

PACS 95.30.Qd • 98.38.Ly • 95.10.Bt • 97.10.Me • 
97.60.Bw 



1 Introduction 

The self-similar gas dynamics in spherical symmetry in- 
volving self-gravity and thermal gas pressure has been 
studied over past several decades with complementary 
perspectives and various applications. In astrophysical 
and cosmological contexts, Larson (1969) and Penston 
(1969a, b) independently studied self-similar flow solu- 
tions in a self-gravitating gas. For modelling star for- 
mation processes in a molecular cloud and in contrast 
to the earlier results (Larson, 1969; Penston, 1969a, b), 
Shu (1977) explored self-similar collapse behaviours of an 
isothermal gas and obtained the central free-fall asymp- 
totic solution and the static solution at large radii ini- 
tially. Shu (1977) focused on the expansion-wave collapse 
solution (EWCS) by numerically joining the inner free- 
fall collapse solution with the outer static isothermal 
sphere solution (i.e., the outer part of a singular isother- 
mal sphere). In different notations, Hunter (1977) con- 
structed complete isothermal self-similar solutions cross- 
ing the isothermal sonic critical line once by matching 
solutions in the speed-density phase diagram. The varia- 
tion trend of a spiral pattern in the speed-density phase 
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diagram suggests that there may exist infinitely many 
discrete solutions. Hunter succeeded in connecting the 
two parts {t <0 and t > with t being the time) of flows 
smoothly. In the same model framework, Whitworth and 
Summers (1985) identified and distinguished such sonic 
critical points as nodal and saddle points through a sys- 
tematic analysis of the behaviours in the neighborhood 
of the isothermal sonic critical line, noted the numeri- 
cal stability issue of integration directions in the vicinity 
of the sonic critical line, and suggested two-dimensional 
continua of solutions with weak discontinuities across 
the sonic critical line. Hunter (1986) promptly examined 
their continua of solutions and pointed out the weak dis- 
continuity in their solutions (Whitworth and Summers, 
1985), suggesting that these solutions may be unphysical 
for being unstable (see also Lazarus 1981 for more de- 
tails). Moreover, Hunter (1986) proposed to use as many 
expansion terms as possible for a numerical integration 
away from nodes along the sonic critical line. Lou and 
Shen (2004) emphasized that the EWCS solution (Shu 
1977) being static sufficiently far away only represents a 
special limiting case of a more general class of constant 
speed solutions at large x and constructed global semi- 
complete solutions for envelope expansion with core col- 
lapse (EECC), connecting the inner free-fall asymptotic 
solutions with asymptotic flow solutions at large x us- 
ing the similar matching procedure of Hunter (1977). In 
particular, Lou and Shen (2004) constructed EECC solu- 
tions crossing the isothermal sonic critical line twice with 
radial similarity oscillations in the subsonic regime and 
with the divergent free-fall asymptotic behaviour in the 
limit of small x. We note in passing that steady or self- 
similar accretions of dark matter under self-gravity might 
be relevant in understanding the formation of a few re- 
cently reported supermassive black holes (SMBHs) in the 
early Universe (Hu et al., 2005). Lou (2005) outlined 
a two-fluid similarity dynamics which is fairly similar 
to the isothermal model mentioned above, to model the 
gravitational coupling between a dark matter halo and 
a hot interstellar gas medium. Furthermore, by incorpo- 
rating a random magnetic field into the model analysis 
(Lou, 2005), it is possible to set up a model framework to 
further examine synchrotron radio emissions and mag- 
netic Sunyaev-Zel'dovich effect in galaxy clusters (Hu 
and Lou, 2004). Also, the stability problem has been 
tackled (Ori and Piran, 1988; Hanawa and Matsumoto, 
1999, 2000; Hanawa and Makayama, 1997) for the Hunter 
type isothermal solutions. However, because of technical 
issues, the free-fall solutions (Shu, 1977) remain to be 
analyzed in this aspect. Semelin et al. (2001) took into 
account of viscosity to examine the flow stability prob- 
lem. 

Physically, we imagine that shocks occur naturally 
under various astrophysical flow situations (see, e.g.. Ken- 
nel and Coroniti 1984 for a model of the Crab Nebula 
involving MHD pulsar wind shocks and Bagchi et al. 
2006 for observational evidence of MHD galaxy cluster 



wind shocks). In contexts of star formation, self-similar 

isothermal shock flows have been investigated and ap- 
plied to astrophysical systems such as Bok globules and 
the so-called 'champagne flows' in Hjj regions (Shu et 
al., 2002; Tsai and Hsu, 1995; Shen and Lou, 2004; Bian 
and Lou, 2005). We briefly touch upon the subject of 
self-similar MHD shocks here, because our MHD results 
(Yu and Lou, 2005) are candidate solutions for both the 
upstream and downstream regions across an MHD shock. 
The self-similar polytropic MHD shock flows will be pre- 
sented elsewhere in more details (see, e.g., Yu et al., 2006 
for self-similar isothermal MHD shock flows and Lou and 
Wang, 2007). 

In polytropic MHD collapses and outflows, the poly- 
tropic index 7 represents a gross simpliflcation from com- 
plicated physical processes including possible nuclear re- 
actions, energy transport, neutrino transport and elec- 
tron capture and so forth in various contexts (see, e.g.. 
Bouquet et al., 1985 and Yahil, 1983). In parallel with 
earlier results under the isothermal approximation, the 
polytropic treatment has also been pursued in various 
astrophysical contexts including supernovae. Goldreich 
and Weber (1980) studied the 7 = 4/3 case with a focus 
on the homologous collapse of the inner core of a pro- 
genitor star prior to the emergence of a rebound shock 
(e.g., Lou and Wang, 2006, 2007). Yahfl (1983) noted 
limitations of Goldreich- Weber model, studied the poly- 
tropic gas dynamics within the polytropic index range of 
6/5 < 7 < 4/3, and discussed possible applications to the 
pre-catastrophe as well as post-catastrophe phases sepa- 
rately. A few years earlier, Cheng (1978) investigated the 
polytropic hydrodynamics with 1 < 7 < 5/3, focussing 
on the initial distribution of mass density p with a radial 
scaling of p ^ r". For a polytropic gas, Cheng (1978) 
derived the inner free-fall solution and generalized the 
isothermal EWCS to the polytropic EWCS. Bouquet et 
al. (1985) introduced a dual space of parameters as well 
as the systems I and II for the polytropic gas dynamics 
with 1 < 7 < 4/3, and analyzed properties of the sonic 
critical curve in details. Suto and Silk (1988) performed a 
similarity transformation for the polytropic gas dynam- 
ics and, from the resulting nonlinear ordinary differential 
equations (ODEs), obtained regular and asymptotic so- 
lutions as well as numerical solutions with n = I and 
"f — 2 ~ n for 1 < 7 < 4/3. By taking into account of 
the total radiative emissivity from the gas in the form 
of e ~ /9°r'*+^, Boily and Lynden-BeU (1995) replaced 
the polytropic equation of state with an energy equa- 
tion and discussed physical, mathematical and numeri- 
cal properties of such radiative self-similar gas dynam- 
ics. McLaughlin and Pudritz (1997) considered the lim- 
iting case of the so-called logotropic gas dynamics also 
involving sonic critical points and obtained expansion- 
wave collapse solutions in their analysis. In the context 
of star formation, Fatuzzo et al. (2004) studied a com- 
bination or a transition of an initial 'equation of state' 
P oc and a later dynamic equation of state P oc p'^, 
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and noted that in the special case oi F — 1 , asymptotic 
solutions may have constant flow speed at far away re- 
gions, analogous to the earlier results of Lou and Shen 
(2004) and Shen and Lou (2004). Fatuzzo et al. (2004) 
carried out their analysis for the case of 7 < 1 without 
involving the sonic critical curve. 

Magnetic fleld can be extremely important in many 
astrophysical processes on different scales and in partic- 
ular, for star formation activities at various stages (e.g., 
Shu et al., 1987). In mostly neutral gas medium, such as 
molecular clouds and cores etc., the magnetic field can 
only couple to the gas medium if MHD wave frequencies 
are less than the ion-neutral collision frequency, corre- 
sponding to a lower limit on the wavelength that MHD 
waves need in order to propagate in such a magnetized 
cloud (e.g., Myers, 1998). At the end of stellar evolution, 
magnetic fields are observed to exist in various stellar 
systems such as the well-known Crab Nebula (e.g., Wolt- 
jer, 1957, 1958a, b; Kennel and Coroniti, 1984a, b; Wil- 
son et al., 1985; Lou, 1993; Wolf et al., 2003). Chiueh and 
Chou (1994) discussed the gravitational collapse of an 
isothermal magnetized gas cloud, including the magnetic 
pressure force term in the radial momentum equation to- 
gether with the magnetic induction equation. They as- 
sumed a randomly distributed magnetic field such that a 
quasi-spherical symmetry is sustained during the MHD 
similarity evolution of a gas cloud on large scales. Mag- 
netic tension force was ignored in their formulation j^Yu 
and Lou (2005) approached the same physical problem 
yet with a different formulation and discussed MHD con- 
sequences of a random magnetic field. Using the frozen- in 
condition on magnetic field, Yu and Lou (2005) managed 
to reduce the three apparently coupled nonlinear MHD 
equations to two key coupled nonlinear MHD ODEs in 
an equivalent manner. Yu et al. (2006) further explored 
various self-similar isothermal MHD shocks under the 
approximation of large-scale quasi-spherical symmetry. 

Self-similar gas dynamics for stellar collapse prob- 
lems under self-gravity and thermal pressure have been 
studied extensively from various perspectives. In this pa- 
per, we construct similarity MHD solutions to explore 
nonlinear effects of a random magnetic field. In general, 
MHD similarity collapses and outflows evolve nonlinearly 
by simply changing certain profile scalings of the en- 
closed mass, gas mass density, radial flow speed, and 
mean transverse magnetic field energy density without 
changing their shapes. Hydrodynamic simulations have 
pointed to possible similarity evolutions as all sorts of 
transients peter out in time (see, e.g., Bodenheimer and 
Sweigart, 1968 and Foster and Chevalier, 1993). The 
well-known example is the Sedov- Taylor similarity blast 
waves resulting from a point explosion (Sedov, 1959; Lan- 
dau and Lifshitz, 1959; Barenblatt and Zel'dovich, 1972). 



The magnetic induction equation in Chiueh and Chou 
(1994) involves typos, both incorrect in its original form and 
inconsistent with their later nonlinear ODEs after the simi- 
larity MHD transformation. 



Self-similar fiow solutions have been explored in different 
geometries (Fillmore and Goldreich, 1984; Hennebelle, 
2003; Terebey et al, 1984; Inutsuka and Miyama, 1992; 
Shadmehri, 2005; Krasnopolsky and Konigl, 2002; Shen 
and Lou, 2006) and we here work with the quasi-spherical 
geometry mainly by a more phenomenological consider- 
ation. For the example of the Crab Nebula projected 
onto the plane of the sky, a quasi-spherical morphology 
(more or less elliptical in reality) has been sustained on 
large scales. Another example is the Cassiopeia A super- 
nova remnant (resulting from a type II supernova explo- 
sion presumably) which, projected onto the plane of sky, 
appears more or less round with a central neutron star 
manifested as a bright X-ray point. There are also ex- 
amples of more or less round planetary nebula systems 
where magnetic fields, be they weak or strong, are also 
likely involved. We invoke these morphological examples 
involving magnetic fields to justify an MHD collapse and 
expansion problem with a quasi-spherical symmetry on 
large scales as a first approximation. We also assume that 
small-scale deviations from the quasi-spherical symmetry 
is relatively insignificant in large-scale MHD, i.e., small- 
scale transverse fiow components are random and may be 
neglected to simplify the mathematical treatment. Since 
magnetic field strengths can be significant in various as- 
trophysical systems (Yu and Lou, 2005), we should take 
into account of the MHD infiuence in the evolution of a 
magnetized gas cloud or a magnetized star (Lou, 1993, 
1994) as well as MHD gas systems on much larger scales. 

For a random magnetic field in a cloud, we envision a 
simple 'ball of thread' scenario in a vast spatial volume 
of gas medium. A magnetic field line follows the 'thread' 
meandering within a thin spherical 'layer' in space in a 
random manner. In the strict sense, there is always a 
random weak radial magnetic field component such that 
random magnetic field lines in adjacent 'layers' are ac- 
tually connected throughout in space. By taking a large- 
scale ensemble average of such a magnetized gas system, 
we are then left with 'layers' of random magnetic field 
components transverse to the radial direction. Having 
gone thus far in our idealization, we would admit that 
the magnetic fields in the Crab Nebula, the SNR Cas A 
as well as several round planetary nebulae may not be 
fully represented by of "ball of thread" scenario. What 
we have been trying to emphasize is the large-scale quasi- 
spherical geometry of magnetized astrophysical systems 
rather than detailed magnetic field configurations. We 
note also that, in our model, the "ball of thread" sce- 
nario is mainly for the transverse magnetic field effect 
on average, while the MHD effect of a weak radial mag- 
netic field may be negligible. As a matter of fact, we 
will still need further observational information to infer 
whether our magnetic field configuration can roughly de- 
scribe some round-shaped morphologies of astrophysical 
systems. 

In reference to the recent isothermal self-similar MHD 
analysis (Yu and Lou, 2005), we show in this paper that 
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an isothermal similarity MHD treatment can be natu- 
rally extended to a magnetized polytropic gas in a sys- 
tematic manner. Parallel to the self-similar transforma- 
tion for relevant variables (Suto and Silk, 1988) with 
an additional transformation for the transverse magnetic 
field, we derive three apparently coupled nonlinear MHD 
ODEs, as in the case of an isothermal magnetofluid (Chi- 
ueh and Chou, 1994). The major technical difference in 
our polytropic MHD formalism is that these three ODEs 
can be readily reduced to two key ODEs of MHD by in- 
voking the frozen-in condition on magnetic field (Yu et 
al., 2006; Lou & Wang, 2007). This frozen-in conditioi|^ 
(21) naturally leads to an integration constant h denot- 



ing physically the ratio of the magnetic energy density 
to the self-gravitational energy density and significantly 
reduces the complexity of analyzing the nonlinear simi- 
larity MHD problem. Although only a change of equation 
of state is made in our current formulation as compared 
to the isothermal treatment (Yu and Lou, 2005), several 
distinct differences arise. For example, the magnetosonic 
critical curve now shows qualitatively different asymp- 
totic behaviours as compared to an isothermal gas, both 
in unmagnetized and magnetized cases (Lou and Wang, 
2006, 2007). By increasing the polytropic index 7, we find 
qualitative differences in reference to the case of a smaller 
7. Most importantly, we found a novel asymptotic non- 
linear MHD solution near the central core or at later 
time and constructed semi-complete MHD similarity so- 
lutions using this asymptotic solution. We have also dis- 
covered the so-called 'quasi-static' asymptotic polytropic 
MHD solution behaviours (see Lou and Wang, 2006 for 
polytropic hydrodynamic asymptotic solutions) . We here 
focus on the MHD case and provide a description in Ap- 
pendix|G] A more detailed analysis and astrophysical ap- 
plications of this MHD asymptotic solution can be found 
in Lou and Wang (2007). 

Motivated by potentially wide astrophysical appli- 
cations, the main purpose of this paper is to present 
possible similarity solutions from the nonlinear MHD 
ODEs, distinguish the asymptotic behaviours of differ- 
ent types including the eigensolutions across the magne- 
tosonic critical curve, and construct global semi-complete 
solutions numerically. Our analyses and results here serve 
as the theoretical basis for further specific astrophysical 
MHD applications. We provide the background informa- 
tion in Section [1] as an introduction. Section [2] contains 
the basic MHD formulation of the problem and section 
[3] presents the mathematical analysis. Section |4] mainly 
describes numerical results, including the magnetosonic 
critical curves, similarity MHD solutions without cross- 
ing the magnetosonic critical curve, similarity solutions 
crossing the magnetosonic critical curve once and twice. 
In both analytical and numerical analyses, we focus on 
differences between the cases with or without magnetic 

^ By combining conservations of mass and magnetic flux, 
we can readily derive equation (21 1 or in dimensional form 
< Bt > l{p^r^) =consant. 



field and between weak and strong magnetic field. We 
also compare the case in which 7 is almost unity to that 
in which 7 is larger than one, and further discuss differ- 
ences between a nearly isothermal polytropic gas and an 
exact isothermal case. 



2 Similarity MHD Flows 

In this section, the basic MHD formulation of the sim- 
ilarity problem is presented and the approximation of 
quasi-spherical symmetry is discussed (Appendix [A| . 



2.1 MHD Formulation of the Problem 

Under the assumptions of a random magnetic field on 
smaller scales, the approximation of quasi-spherical sym- 
metry and the ideal MHD treatment, the dynamics of a 
polytropic magnetized gas in spherical polar coordinates 
(r, 0, 0) is described by the following equations: 



dt 
dM 



)_d_ 

j,2 Qj, 



(r^pu) = 



dM ^ 
at or 



dM 
dr 

fdu 



Airr p 



du 
dr 



(1) 
(2) 
(3) 



GMp 



p = Kp^ 



dp 
dr 
_d_ 

dr \ 8tt J Airr 

2r'<Bf>^ = 0, 
dr 



(4) 



(5) 



(6) 



where G ~ 6.67 x 10~^ g~^ cm~^ s~^ is the gravitational 
constant, p{r, t) is the gas mass density, M(r, t) is the en- 
closed gas mass within radius r at time t, u{r, t) is the 
bulk radial fiow speed, and < > is the mean square of 
the random transverse magnetic field Bj proportional to 
the magnetic energy density associated with the random 
transverse magnetic field. In the conventional polytropic 
equation of statc[^ ^ , the coefficient k remains constant 
globally, independent of both r and t. The case of 7 = 1 is 

^ The condition of specific entropy conservation along 
streamlines would be more general and will be considered 
in a separate paper. 
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only a special case corresponding to an isothermal mag- 
netized gas (Yu and Lou, 2005; Yu et al., 2006). The 
Poisson equation relating the mass density and the grav- 
itational potential is automatically satisfied under the 
quasi-spherical symmetry. In the above equations, the ra- 
dial momentum equation Q involves the magnetic pres- 
sure and tension forces on the right-hand side (RHS), 
and equation (|5| is derived from the magnetic induction 
equation along with certain simplifications (Yu and Lou, 
2005); these two equations will be further discussed and 
analyzed in the next subsection. Compared to the work 
of Chiueh and Chou (1994), the formulation here is dif- 
ferent by keeping the magnetic tension force and by deal- 
ing with a conventional polytropic gas. For the problem 
outlined above, the reader may consult relevant refer- 
ences (Shu 1977; Suto and Silk, 1988; Chiueh and Chou, 
1994; Lou and Shen, 2004; Bian and Lou, 2005; Yu and 
Lou, 2005; Yu et al., 2006; Lou & Wang, 2006, 2007). We 
focus on the semi-complete solution space < t < -l-oo 
rather than the complete solution space as introduced by 
Hunter (1977)0 

For a conventional polytropic gas with 7 = 2 — n 
where 7 and n are defined by equations (15 1 and ([9]), 



respectively (see subsection 3.1 ), a combination of equa- 
tions ([1]) — (|5| together with equation (15) leads to the 
MHD energy conservation equation as 



d_ 
di 



pu 



p 



(7-1) 

]_d_ 
1 d ( 



(9^> 

dr 



<Bi > 



{8n) 



r^up 



IP 



r u 



(7-l)p 



(7) 



where —d^P/dr — —GMjr^ and ^ is the gravitational 
potential (Fan and Lou, 1999). This MHD energy conser- 
vation equation reduces to the isothermal cases (Lou and 
Shen, 2004; Yu and Lou, 2005) by taking the L'Hopital 
rule with respect to 7 in the limit of 7 ^ 1. The mag- 
netic energy density and Poynting flux density associated 
with < B\ > can be readily identified in the MHD en- 
ergy conservation equation ([7| . With the quasi-spherical 
symmetry, the divergence term containing (p in MHD en- 
ergy conservation equation ([7| vanishes (Lou and Shen, 
2004). 



2.2 Comments on the MHD Formalism 

We here briefly comment on the basic MHD formulation, 
the quasi-spherical symmetry, the three-dimensional ran- 

* By the time-reversal invariance, the correspondence be- 
tween a complete solution (Hunter, 1977) and semi-complete 
solutions has been shown explicitly in Lou and Shen (2004) 
by concrete examples. 



dom flow fluctuations on small scales and the physical 
basis for the magnetic force density and the magnetic 
induction equation. 

Our MHD model describes a self-gravitating gas cloud 
embedded with a magnetic field presumed to be random 
and tangled on small scales. On large scales, the gas mass 
density, the gas thermal temperature, the thermal pres- 
sure and the entropy are all taken to be quasi-spherically 
symmetric. The magnetic field distribution is presumed 
to be completely random in space such that in a small 
volume (an infinitesimal volume r^drd^d^) the magnetic 
field is effectively represented by the mean square aver- 
ages of < > and < Bf >, proportional to the radial 
and transverse magnetic energy densities respectively. In 
terms of the magnetic pressure and tension forces for the 
large-scale MHD, < Bf > plays the dominant role on 
the dynamics of a magnetized gas cloud as compared to 
< i?^ >. As the magnetic field is randomly distributed 
with a quasi-spherical symmetry, the bulk gas flow ve- 
locity remains grossly spherically symmetric and can be 
characterized by the bulk mean radial flow speed Vr] the 
transverse component of the flow velocity vg and are 
relatively small and may be neglected in the first approx- 
imation. These transverse components should be part 
of Alfvenic fluctuations corresponding to magnetic field 
fluctuations about the mean configuration; thus the more 
random the magnetic field fluctuations are, the better the 
approximation becomes. 

The physical concept of a quasi-spherical symmetry 
for a magnetized gas cloud or a magnetized progeni- 
tor star is only valid for MHD processes of sufficiently 
large scales. Here, 'large scales' are obviously in contrast 
to 'small scales' on which magnetic fields are presumed 
completely random locally in our MHD model frame- 
work. In the strict sense, an exact spherical symmetry 
is impossible due to the very nature of a magnetic field. 
However, a quasi-spherical symmetry may be sustained 
for large-scale MHD processes. We invoke the projected 
quasi-elliptical shape of the Crab Nebula and the pro- 
jected more or less round remnant of the Cassiopeia A 
supernova as empirical supports for this notion of quasi- 
spherical symmetry. Although it is not yet obvious that 
the actual magnetic field can be largely approximated 
by our 'ball of thread' scenario, from the morphology of 
these systems we suggest a globally random magnetic 
field distribution as a plausible yet tractable starting 
point. In this scenario, small-scale random flow velocities 
are ignored as compared to systematic radial flows. Qual- 
itatively speaking, this perspective is justifiable when a 
random magnetic field is weak. In our model analysis, 
we sometimes do encounter situations of strong magnetic 
fields especially for accretions towards a central compact 
object. In such a case, one should really view our asymp- 
totic MHD solutions as indicating a gross trend of varia- 
tion that is bound to be destroyed by non-spherical and 
transient MHD processes sufficiently close to the central 
compact object. Upon impacting onto a central compact 
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object, we expect the emergence of a strong radiating 
MHD shock travehng outward slowly in a self-similar 
manner (Shen and Lou, 2004; Yu and Lou, 2005; Yu 
et al., 2006). Practically, we can apply our large-scale 
self-similar MHD solutions of quasi-spherical symmetry 
outside this radiating MHD shock. Within this quasi- 
spherical MHD shock, the core gravity can be strong 
enough to more or less hold on the strongly magnetized 
plasma. By the naive solar analogy, we readily imagine 
that sporadic violent 'flares' or 'coronal mass ejections' 
may erupt from the strongly magnetized central core re- 
gion and can even break into the 'self-similar' and 'quasi- 
spherical' domain of magnetized accreting flows. 

Given local random magnetic fields on small scales 
in a gas medium, three-dimensional MHD flows on small 
scales are naturally expected because of the unbalanced 
magnetic tension force here and there. Except for very 
special self-similar radial flow situations of zero trans- 
verse flows yet with a three-dimensional magnetic fleld 
(see Low 1992), we do generally expect small-scale three- 
dimensional random flow fluctuations associated with 
the large-scale mean radial flow. By our assumption of 
a randomly tangled magnetic fleld, such flow fluctua- 
tions are more or less confined or trapped locally and 
advected by the mean radial MHD fiow on large scales. 
In short, we do not expect mean flows transverse to the 
radial direction on large scales in our scenario. By in- 
tuition, we expect transverse flows caused by the ran- 
dom magnetic tension force, yet such flows will remain 
locally confined due to the local random magnetic field 
and hence be small as compared with the bulk quasi- 
spherical radial fiow speed. As suggested by Zel'dovich 
and Novikov (1971), an isotropic magnetic pressure is 
expected from a completely random magnetic field on 
small scales. We follow this basic concept and also in- 
clude the radial magnetic tension force which is non- 
negligible in our formulation. Physically, one may view 
such small-scale random fiow fiuctuations as turbulence 
and for simplicity, we have ignored the effects of the ef- 
fective turbulent pressure, viscosity and resistivity etc. 
(see, e.g., Lou and Rosner, 1986) in this formulation. 
In our model framework, if the MHD turbulence will 
attribute to random fields and fiuctuations on smaller 
scales, the turbulence scale itself should be small enough 
such that the largest turbulence scale is small compared 
to the overall quasi-spherical geometry. 



Magnetized Polytropic Gas Cloud 

To seek self-similar solutions to the MHD equations, we 
introduce an independent similarity dimensionless vari- 
able X and presume that dependent physical variables 
are given by the following similarity forms accordingljj^ 



r — ax, u = bv, p — ca, p — df3, M — em, < >— fw , 

(8) 

where the six scaling factors a{t) through f{t) are func- 
tions of time t only and are defined by 



6 = /fci/2r-i , 

" (3n-2)G ' 



1 



G 



(9) 



Here k and n are two constant parameters. As func- 
tions of X only, v{x), a{x), (3{x), m{x), and w{x) are 
the reduced forms of radial fiow speed, gas mass density, 
gas pressure, enclosed gas mass, and magnetic energy 
density (associated with the averaged random transverse 
magnetic field) , respectively. With this self-similar MHD 
transformation, equations ^ and (|3| lead to an alge- 
braic expression for 'm(x) in terms of a(x) and v(x), viz. 



m — ax'^{nx — v) , 

and an ODE for a(x) and v{x) 

[nx — v)a ~ av — ~2 a 



(10) 



(11) 



where the prime ' denotes the differentiation with respect 
to X. Relation (10 1 leads to the important inequality 



nx — V > 



(12) 



for a positive gas mass density as noted repeatedly in 
figure displays presently. In our later analyses, this in- 
equality is a key constraint on choosing relevant physical 
solutions. From equation ([5|, one obtains 



, 2v~{A~ 2n)x 
[nx — v)w — Zwv = w 



(13) 



3 Model Analysis 

With the basic ideal MHD model qualified and the for- 
mulation established, we now perform the analytical and 
numerical analyses in order. In the next section, we present 
the results of numerical exploration. 



3.1 Self-Similar MHD Processes in a 



for the reduced dependent variables w{x) and v{x). This 
constraint for the reduced magnetic energy density is 
fairly similar to equation (111. By equation ([3]) , one ob- 
tains the reduced radial momentum equation 



\nx — v)v -I- — = —\n—\)v^ 



2a 



in- 2 



(14) 



In terms of the self-similar transformation, the magnetic 
field term here distinguishes ours from that of Suto and Silk 
(1988). 
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For a generalized polytropic equation of state with a 
polytropic index 7 (Suto and Silk, 1988), we simply have 



P 



(15) 



where k may be time dependent in general. For n — 2— 7, 
we have a constant k = /c(47rG)^~^ and equation (15 1 is 



an equation of state for a conventional polytropic gas. 
A combination of equations (14 1 and (15) leads to 



7a 



7-2, 



a' — (nx — v)v' + 
= — (n — 1)1; — 



w' 
2a 
{nx — v) 



xa {3n — 2) 



(16) 



and equations (11 



(131 and (^16 are the three MHD 



similarity ODEs describing polytropic magnetized gas 
flows with a quasi-spherical symmetry. The three nonlin- 
ear MHD ODEs are similar to those of Chiueh and Chou 
(1994) with the key differences in the adopted equation of 
state and in keeping the magnetic tension force term (see 
equations [T| to |5| . By taking relevant limits as necessary 
checks, these equations are consistent with those of Shu 
(1977), Suto and Silk (1988), Lou and Shen (2004), Yu 
and Lou (2005), Yu et al. (2006), Lou and Wang (2006, 
2007) as expected. 

The Alfven speed in this formulation is defined by 



_ ( <Bf > 
Airp 



1/2 



(17) 



and the sound speed s in the polytropic gas determined 
by the equation of state with 7 = 2 — 71 (i.e., the poly- 
tropic state equation in the usual sense) is simply 



dp 
dp 



1/2 



(18) 



Thus the ratio of the Alfven wave speed to the polytropic 
sound speed becomes 



VA _ / w \ 
s y^a^ J 



1/2 



(19) 



consistent with the isothermal 7=1 case (Yu and Lou, 
2005; Yu et al., 2006). 



3.2 Reduction of Nonlinear MHD ODEs 

One can readily reduce the three coupled nonlinear MHD 



ODEs to two. From equation (111, one obtains 

v' — 2{x — v)/x 



a' 



a 



nx — V 



and from equation (13 1, one gets 
w' 2v' + 2[2v - {'i-2n)x\/x 



The above two equations lead to a differential relation 

^ = ^ + ^ (20) 
wax 

which immediately gives a simple integral of 

(21) 



I, 2 2 
w — ha X 



where h is an integration constant providing a measure 
for the magnetic field strength. Integral (21 1 gives a new 



parameter /i of a magnetized gas cloud besides n and 7, 
and reduces the three coupled nonlinear MHD ODEs to 
two. Expressed explicitly in physical quantities, h is 



h = 



representing the ratio of the magnetic energy density to 
the self-gravitational energy density. This simplification 
reduces tremendously complications in numerical MHD 
exploration and physically represents the frozen-in con- 
dition on magnetic field (Yu and Lou, 2005; Yu et al., 
2006) . The procedure of numerically constructing global 
MHD solutions and matching the solutions across the 
magnetosonic critical curve can be carried out similar 
to that of Lou and Shen (2004) . This reduction of three 
coupled nonlinear MHD ODEs to two is parallel to the 
isothermal case (Yu and Lou, 2005; Yu et al., 2006). 
Substitution of equation ( 21 ) into equation ( 16 1 gives 

(70'''^^ + hx'^)a' — [nx — v)v' 



-(n — l)u - 



2hx 



'in -2 



(22) 



Equations (11) and (^22^ together lead to the two coupled 



nonlinear MHD ODEs in the forms of 
2 



a = a 
2{^x — v)(nx — v) 



(n-l)i;+ 2hx+ \a 



3n- 2 

a{nx — vf' — "fa'^ — ha^x"^ 

I {tlx 

v' = l{n- I) [av{nx - v) + 2ha^x^] + ^ ^a~ 



(23) 



(3n - 2) 



-270;' 



a{nx — u)'^ — ja"^ — hcP'x^ 



(24) 



We now introduce simplifying notations as follows 



A(q!, x) = 



(n — l)u 



nx ~ V \ (x — v)(nx — v) 

2/ixH a -2^ ' 

in -2 



X 



2^21 



V{a, V, x) = {n ~ 1) \av{nx — v) + 2ha x 



{nx — v)' 
(3n - 2) 



a ~ 2ja' 



,{x~v) 



(25) 



(26) 
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and 



X{a, V, x) = a{nx — v)"^ — ^a'^ 



to transform equations ( 23 1 and ( 24 1 into 
da; 



da 



du 



Xia^v^x) A(a,ti,x) l/(a,ti,x) 
where ^ is a new dependent variable. 



(27) 



(28) 



The two coupled nonlinear MHD ODEs ( 23 1 and ( 24 1 



are analyzed to determine the magnetosonic critical curve 
and asymptotic solution behaviours near the magnetosonic 
critical curve, and can be integrated numerically using 
the standard fourth-order Runge-Kutta scheme (e.g., Press 
et al., 1986). Along with the similarity MHD transforma- 
tion, these two coupled nonlinear MHD ODEs describe 
an important subset of MHD solutions to the original 
nonlinear partial differential MHD equations ([T]) — ([5|. 

Using the above simplification, the ratio of Alfven 
wave speed va to the gas sound speed s becomes 



vaI s = (ha X /7) 
for 7 + n = 2. 



a/2 



(29) 



3.3 Singular Surface and Magnetosonic Critical Curve 



Given parameters n, 7 and h in MHD ODEs (23 1 and 
(24), there exists a characteristic surface in the (iJ^ a, a;) 



space on which the denominators on the RHSs of both 
equations vanish (Whitworth and Summers, 1985). Phys- 
ically, we have averaged over small-scale MHD fluctua- 
tions in our formulation. Therefore, the magnetosonic 
critical point or curve should correspond to a layer of 
a thickness comparable the mean scale of MHD fluctu- 
ations. In other words, our model analysis relates to an 
averaged condition in the actual MHD flow. Mathemati- 
cally, this singular surface is determined by computing v 
from specific x and a in the ranges of < a; < +00 and 
< a < -|-cx), namely. 



(30) 



in which one should pick up the upper minus sign in or- 
der to satisfy the physical constraint of na; > ?; and thus 
to ensure mix) > 0. The solutions of the two coupled 
ODEs (23 1 and (24) cannot cross the singular surface 



unless they cross it at points along the so-called mag- 
netosonic critical curve where both the numerators and 
denominators vanish simultaneously. Along this magne- 
tosonic critical curve, the derivatives of the dependent 
variables v{x) and a{x) can be calculated from equations 
(23) and (24) using the L'Hopital rule. Mathematically, 



this magnetosonic critical curve is defined by the follow- 
ing pair of equations 



{nx — w)^ — ja^ ^ + hax"^ 



(31) 



nx ^ v\ 2{x — v){nx — v) 



{n-l)v+{2hx+ 

\ in — 1 ) X 

these two equations immediately give 

V = nx^ (70;'''""'" -I- hax^Yl"^ , 



0; (32) 



± 



1 



(3n - 2) 



27 



(7a''' 
7-1 



+ /iax2)i/2 
n(n — l)a: , 



(33) 



with the latter leading to a quadratic equation of x^ as 
Axx^ + Bxx^ + Ci = , (34) 
where the coefficients A\, B\ and C\ are deflned by 



n-\ 



Bi = 7a('^~i) 



(3n - 2) 



ri — 1 



ha — r?(n — 1)^ 



(3n- 2) 



-f 4n(n- 1) 



(35) 



If one substitutes equations ( 32 1 and ( 33 1 into equa- 
tion (24), the numerator vanishes. By the physical con- 
straint of m(a:) > 0, the lower minus sign in equation 
(33 1 will be ignored, even though mathematically, it may 



represent a new branch of the critical curve should this 
branch do exist. The above expressions of the magne- 
tosonic critical curve appear far more complicated than 
the isothermal case (Shu, 1977; Lou and Shen, 2004; Yu 
and Lou, 2005; Yu et al., 2006; Lou and Gao, 2006), as 
a result of a polytropic gas under the influence of a ran- 
dom magnetic field characterized by a constant h and 
a reduced magnetic energy density wi^x). In reference 
to earlier results of determining ?;(a;) and a(x) in terms 
of X along the critical curve, the most straightforward 
procedure one can take in the current polytropic MHD 
problem is to first determine x from a given a and then 
obtain the corresponding v. This is somewhat unusual in 
determining the magnetosonic critical curve for a given 
sequence of a values. The additional constraints for the 
magnetosonic critical curve in the semi-complete space 
are a; > and a > besides equation ( 34 1 . Physically, we 
are interested in the parameter regime of na; — f > such 
that m,(x) > 0. It is obvious that we always have Ci < 
in definition ( 35 1 . For different a values and depending 
on the values of coefficients Ai, Bi and Ci, there are 
three possible cases listed below. 

Case I: subcase (i) of both Ai < and _Bi < or sub- 
case (ii) a negative determinant Bf — AAiCi < 0; there 
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is then no positive root for satisfying quadratic equa- 
tion ( 34 1 and therefore there is no point along the mag- 



netosonic critical curve corresponding to such a range of 
a values. 

Case II: with Ai < 0, Bi > and a non-negative 
determinant — AAiCi > 0, there are two positive 



roots for x satisfying quadratic equation ( 34 1 ; there are 



thus two points on the magnetosonic critical curve cor- 
responding to such a range of a values. 

Case III: with Ai > 0, there is only one positive root 



for X satisfying quadratic equation ( 34 1 and there is 



thus one point on the magnetosonic critical curve corre- 
sponding to such a range of a values. 

Once an x value is determined for a given a > by 
the above procedure, we readily obtain the corresponding 
V value. According to equation (33 1 for 



27a''' ^/x — n(n^l)x 
[n-l + a/(3n-2)] 



> 



(36) 



we should obviously pick up the upper sign (i.e., the 
'+' sign in the second relation and '— ' sign in the first 
relation) in equation (33 1; otherwise, we should pick up 



the other sign accordingly. As the physical constraint 
m{x) > requires that nx — v > 0, inequality (36) sets 
the criterion for a physical solution. 

This sequence of determining the magnetosonic criti- 
cal curve appears more involved than those in the isother- 
mal case without a random magnetic field (Shu, 1977; 
Lou and Shen, 2004). 



3.4 Asymptotic and Global Similarity Solutions 

In order to specify initial or boundary conditions for nu- 
merical integrations, one needs to derive asymptotic sim- 
ilarity MHD solutions. These asymptotic MHD solutions 
also carry their physical implications. It is also possible 
to derive some regular solutions from the coupled non- 
linear MHD ODEs for physical interpretations and for 
reference of numerical results. In general, we have found 
the MHD counterparts of the isothermal asymptotic so- 
lutions (Lou and Shen, 2004), and, in particular, we have 
derived a novel MHD asymptotic solution as well as a 
new class of asymptotic behaviours. We also examine 
the corresponding ratio of the Alfven wave speed to the 
sound speed and the dominant force among the gravity 
force, the thermal pressure gradient force, the magnetic 
pressure gradient force and the magnetic tension force 
for each MHD asymptotic solution. We also analyze the 
behaviour of crossing the magnetosonic critical curve for 
MHD solutions. 

3.4-1 Asymptotic Solutions in the Limit of Large x 



As X approaches infinity in equations ( 23 1 and ( 24 1 and 



with finite a{x) and v{x), one obtains the following pair 



of MHD ODEs to leading orders of large x 

, 2a 

a = , 

nx 



(71 — l)z; 



2h{n- 1) 



27Q! 



7-1 



(37) 



(38) 



371 — 2 

We solve these two asymptotic MHD ODEs to obtain 



V = Bqx^-^/'' 



2h{n- 1) 



3?T. — 2 n 

(2-27-n)/n 



n[2{n -I- 7) - 3J 

where Aq and Bq are two constants of integration. With 
similarity MHD transformation ([s]) and ([9|, this asymp- 
totic MHD solution becomes 



p = ca — cAqx 



2/n _ ^ 



47rG 



(41) 



(42) 



to the leading order of large x, indicating that the gas 
mass density and radial fiow speed profiles are both in- 
dependent of time t at large r. For n = 1, expression 
( 42 1 represents a constant radial fiow speed at very large 



r as emphasized earlier (Shen and Lou, 2004; Lou and 
Shen, 2004) in an isothermal gas. For a{x) and v{x) to be 
non-increasing at large x, this type of asymptotic MHD 
solutions requires that 

max[0, 2(1 - 7)] < n < 1 . (43) 
Note that for h greater than the critical value he, i.e.. 



h > hr — 



2(l-n)(3n-2) 



(44) 



the coefficient of x*-" term in equation (40 1 becomes 



positive, while for h < he this coefficient is negative. The 
presence of this critical value he for ft, is a consequence of 
the polytropic and magnetic nature of our MHD prob- 
lem, which will be further discussed. The corresponding 
reduced mean magnetic energy density w at large x is 



w = Alhx^-^/" 



(45) 



which goes to zero as x ^ +00. With similarity MHD 
transformations ([s]) and we have 



< B^ >= fw = 



^0" . 2/r: 



.2-4/ri 



(46) 



in dimensional form, indicating that the magnetic field 
does not change with time t at large r in a magnetized 
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gas cloud. For a usual polytropic gas with 7 + n = 2, the 
corresponding ratio of the Alfven wave speed va to the 
sound speed s for large x remains constant 



VA 
S 



A2--1 
^0 



1/2 



(47) 



at large x. Ks x appr oache s in finity, the denominators 
of both equations (pSl and (24 1 approach r?ax^ for this 



series of solutions and these solutions will not encounter 
the singular surface w, x) = 0. In the regime of large 
X and for 7 = 2 — n, the dominant forces are both mag- 
netic pressure and gas pressure, and the magnetic pres- 
sure force is stronger than the magnetic tension force 
in magnitude with a ratio of (2/n) — 1. The magnetic 
pressure gradient force points radially outward. 



Since we take 7 > 1, this speed ratio approaches zero 
as a; — > 0+. These asymptotic similarity MHD solutions 
are the free-fall solutions entrained with a random mag- 
netic field. The x— dependence in this limiting behaviour 
is related to the value of n. As x a ppro ache s zero, the de- 



nominators of both equations ( |23| and ( 24 ) approach av^ 
for this series of asymptotic MHD solutions and these 
MHD solutions will not reach the singular surface. Under 
the assumptions above, one further obtains higher-order 
terms for the asymptotic similarity MHD solutions 



a(x) 



(3n- 2)m(0) 
2a?3 

(3ri — 2) 7a: 
(7-1) 2m(0) 



1/2 



(3n - 2)to(0) 



2a;3 



7/2 



(55) 



3.^.2 Asymptotic Solutions in the Limit of Small x 

With the assumptions of ^ a^^^ + hax^ and of 
Q,7-2 ^ XV as X approaches zero in MHD ODEs ( 23 ) and 



(24), we derive the following pair of asymptotic MHD 



equations to leading orders of small x, viz. 

da 2a o? 

Ax X (3n — 2)u ' 

dw a 

da? ~ (3n - 2) ' 

which, by a direct integration, lead to two integrals 
1/2 



(48) 



(49) 



aix) 



(3n-2)m(0) 
2x3 



v(x) = — 



2to(0) 



-,1/2 



(3n- 2)x 



(50) 



(51) 



where m(0) is an integration constant representing the 
core mass at the centre. For this family of asymptotic 
solutions at small x, both assumptions stated at the be- 
ginning of this subsection are satisfied when 



5 

^<3 



and 



2 

"^3 



(52) 



The corresponding reduced magnetic energy density w 



(3n - 2)/(,to(0) 
2x 



(53) 



which diverges as x ^ 0"*" . For a conventional polytropic 
gas with 7 = 2 — n, the corresponding ratio of the Alfven 
wave speed to the sound speed becomes 



VA_ 

s 



(3n- 2)m(0) 



(2-7)/2 



.(37/2-1) 



1/2 



(54) 



v{x) 



2m(0) 
(3ri- 2)a; 

1 72;^ 
(7-l)"^(0) 



1/2 



(3n-2)m(0) 
2a?3 



7/2 



(56) 



with the corresponding reduced magnetic energy density 



, , (3n-2)hm{0) , 

w(x) = — + h-f 

2x 



{3n ~ 2) 



(3n-2)x3 
2to(0) 



(in - 2)m(0) 
2^^ 

1/2 



7/2 



(57) 



In comparison with the isothermal analysis (Whitworth 
and Summers, 1985), this solution appears somewhat dif- 
ferent. This is mainly because Whitworth and Summers 
(1985) considered the isothermal case of 7 = 1; in that 
case the denominators of the second-order terms vanish, 
and one should take the L'Hopital rule with respect to 7 
in order to derive the proper form of expansion solution, 
that is, an a; In a; form for the second-order term. For this 
central free-fall asymptotic solution, the leading force is 
the gravity force, and the magnetic tension force is twice 
the magnetic pressure force in magnitude. The magnetic 
pressure gradient force points radially outward. 

3.^.3 Novel Magnetic Solutions in the Limit of Small x 

For a magnetized gas cloud, it is possible to derive an- 
other MHD asymptotic series solution in which v — 0{x) 
and a — > +00 as x — > 0"*". With the assumption of 



0{x) and a'^ <C a x in equation (24), we have 



v' approaching a constant and obtain 

V — —Cx , 



(58) 



where C is an integration constant. Substituting expres- 
sion (|58|) into equations ( 23 1 and ( 24 ) , we obtain 



[2h + {n + C)/{3n-2)]( 



hx 



(59) 
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Magnetosonic Critical Lines for Different h Values 



Magnetosonic Critical Lines for Different h Values 



7 = 1.01 n = 0.99 



h = 100 



h = W 
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Fig. 1 The magnetosonic critical curves for different h val- 
ues with given parameters n = 0.99 and 7 = 1.01 (7 = 2 — n 
is imposed for a usual polytropic gas). The straight dotted 
line passing through the origin represents nx — v — 0; physi- 
cal MHD flow solutions with positive mass should be to the 
upper right of this straight line. As h increases, the average 
slope d{—v)/dx of an individual magnetosonic critical curve 
increases from negative to positive. The two segments for the 
case of h — with two different line types (i.e., light solid and 
dashed lines) indicate that this curve consists of tw o po rtions 
corresponding to two different roots of equation ( |34[ ). Too 
compact to be seen here, the diverging behaviours ait small x 
for all the magnetosonic curves are shown in Fig. [3] 



Magiretosonic Critical Lines for Different h Values 
7 = 1.25 
n = 0.75 

15 
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h = 10 



h = l 



X 

Fig. 2 The magnetosonic critical curves for different values 
of h with given parameters n = 0.75 and 7 = 1.25 (n = 2 — 7 
for a usual polytropic gas). The light dotted straight line 
nx — V = passes through the origin; the two perpendicular 
light dotted lines represent the abscissa and ordinate axes, 
respectively. Similar to variation trends shown in Fig. [T] as 
h increases (i.e., 0, 1, 10, 100), the average slope d(— i^J/da; 
of an individual magnetosonic critical curve increases from 
negative to positive, and the two segments (i.e., the light solid 
and dashed lines) for the case of h = with two different line 



types correspond to the two sensible roots of equation (341. 
The magnetosonic critical curve with h = 100 (heavy das tied 
lines) has two branches, as indicated in the figure and the 
lower branch beneath nx v — 0. 
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Fig. 3 Enlarged version of Fig. ^ with the same parameters 
yet for much smaller x, intended to clearly show the diverging 
behaviours of the magnetosonic critical curves as a; ^ 0^. 
The horizontal x axis is in a numerical scale of 10""^. 
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Fig. 4 Enlarged version of Figure |2| with the same param- 
eters yet for much smaller x, inteirded to clearly show the 
diverging behaviours of the magnetosonic critical curves in 
the limit of a:: — + 0+. 



-2(n-l)- 



h{3n-2) ■ 
By equation ( [58^ , this in turn requires 
{n + Cf 



2{n-l) + 



h{3n - 2) 



C 



(60) 



(61) 



for consistency. Once this quadratic equation (61 1 for C 



has at least one positive root of C, we obtain at least one 
possible asymptotic MHD solution in the form of 



y{x) = i^^-2-("+C)/[(3n-2)/.] 



(62) 
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v{x) — —Cx , 



(63) 



where D is yet another integration constant. Quadratic 
equation (IM) for C can be readily solved to give 



(3n-2) 



h±{h^-4:hy/^ 



(64) 



and the requirement of nx — u > is satisfied for both 
roots of C. This immediately requires /i > 4 (i.e., a suf- 
ficiently strong magnetic field) for a valid asymptotic 
MHD solution of this kind. In short, this asymptotic so- 
lution for small x is described by 



v{x) = < n — 



(3n - 2) 



h±{h'^- 4h) 



1/2 



- nx-5/2=FV'i^-4/i/(2/i) 



a{x) = Dx 



(65) 



(66) 



and the corresponding reduced enclosed mass m{x) 



(3n-2) ^ 
m(a;) = U 



h±{h^ ~ 4h) 



1/2 



(68) 



The requirement on h is discussed below. For a' ^ a^x^, 
we obtain the following inequality 



7< 



6/i±2(/i2 _4/i)i/2 
5/i± (/i2 -4/i)i/2 



(69) 



For the upper plus signs, if 7 < 6/5, this requirement 
is automatically satisfied, while if 7 > 6/5, this require- 
ment means that h > 4/[l — (57 — 6)^/(2 — 7)^]. For the 
lower minus signs, if 7 < 6/5, this requirement means 
that A<h< 4/[l - (6 - 57)2/(2 - 7)2], while if 7 > 6/5, 
this condition cannot be met, i.e. there does not exist 
such asymptotic solutions with the lower negative signs. 
In the absence of magnetic field with /i = 0, this form 
of asymptotic solution disappears completely. This is a 
brand-new asymptotic MHD solution in a magnetized 
gas cloud, and global semi-complete solutions matching 
this asymptotic solution can be constructed numerically 
(see Figures 15 16 and [Tq] for specific examples). Phys- 



ically, this asymptotic MHD solution describes a much 
compressed accreting nucleus where the magnetic pres- 
sure < > /(Stt) becomes much stronger than the 
thermal gas pressure p to oppose the gravitational col- 
lapse such that the reduced radial inflow speed v{x) ap- 
proaches zero linearly with 0+. Physically, 
we anticipate that a very strong random magnetic field 
confined to a sufficiently small spatial volume would cer- 
tainly destroy the quasi-spherical symmetry at some point 



and drive random flows. We further expect violent and 
sporadic magnetic activities to destroy the similarity evo- 
lution. In spite of all these, we count on the gravity of 
accreted core materials to more or less control a cen- 
tral sphere. In other words, sufficiently far away from 
this central magnetized sphere of influence, we may ig- 
nore feedbacks of central activities and apply our self- 
similar MHD inflow solutions. The scenario envisioned 
here essentially parallels that of a spherical symmetric 
central inflow without magnetic field. Ultimately, there 
must be a central object to confront radial inflows and 
thus destroy the similarity flow evolution. A self-similar 
flow solution is only valid on large scales and outside 
a certain sphere surrounding the core. As x approaches 
zero, the denominators of both equations ( 23 1 and ( 24 1 
for this asymptotic series approach —ha^ 



and these 



MHD solutions do not encounter the magnetosonic sin- 
gular surface. For this asymptotic MHD solution, the 
magnetic pressure, tension and gravity forces are in the 
same order of magnitude, all overpowering the thermal 
gas pressure force, and the magnetic pressure force is the 
strongest. Including one more term in the series expan- 
sion, this novel magnetic asymptotic similarity solution 
appears as 



v{x) 



(3n - 2) 



h±{h'^ -Ah) 



1/2 



_2 (3n -2)[ih-2± 3(/i2 _ 4/i)i/2] 



(6-57)/iT7(/i2-4/i)i/2 

Xa;(2-7) [b/2±VW^/(2h)\ -1 _| 



(70) 



X 6/1 - 4 ± 6(/i2 - 4/i)i/2 4. [5 ± (/j2 _ 4/j)i/2//jj 
x[(6-57)/it(/i^-4/i)^/2] 
X [(6-57)/iT7(/j^-4/i)^/^] 
[(6 - 57)/i ± (2 - -f){h? - 4/^)1/2] 

X2;(7-l)[-5/2Ty)?^/(2/i)]-2 _^ 

3.4^.^ A Singular Global Magnetostatic Solution 



(71) 



For a constant v, equation (11 1 reduces to 



, . , ^{x ~ v) 

(nx — v)a = —2 a 

x 



(72) 



which can be readily integrated for a{x) in the form of 
a{x) = cx-\nx ~ vf-^^'\ (73) 
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Magnetosonic Critical Lines for Different h Values 



b 10 




Fig. 5 Magnetosonic critical curves of a versus x displayed in 
logarithmic scales for different values of h with given 7 — 1.01 
and n = 0.99 as shown in Figure [T] 
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Fig. 6 Magnetosonic critical curves of a versus x displayed in 
logarithmic scales for different values of h with given 7 — 1.25 
and n = 0.75 as shown in Figure [2] 



where c is an integration constant. Along with equation 
( 22 ) , there exists a special singular global magnetostatic 



solution such that v = when rt + 7 = 2, namely 
(2-7)' , (1-7) 



v = 



TO = (2 - 7) 



27(4 - 37) 



7 



27(4-37) 7 
The reduced magnetic energy density w{x) is 



X 



(74) 



w{x) — h 



(2-7)' , (1-7), 



27(4 - 37) 7 



2/(7-2) 



2-4/(2-7) 



(75) 



which diverges as a; ^ 0"*" and vanishes as a; — s- +00. 
Note that for h> he, this global similarity MHD solution 
does not exist. When 7 = 2 — n for a usual poly tropic 
gas, the corresponding ratio of the Alfven wave speed to 
the sound speed becomes a constant 



VA 

s 



1/2 



(2-7)' , (1-7) 



27(4 - 37) 



7 



-1/2 



(76) 



This global similarity MHD solution reduces to that of 
Suto and Silk (1988) for /i = or u> = as expected. This 
is a new singular polytropic magnetostatic solution, con- 
structed in the similar manner as has been done before 
(Shu, 1977; Suto and Silk, 1988; Lou and Shen, 2004). 
When h < he, this solution encounters the magnetosonic 
singular surface at the point 



2n 



2h + n/(3n-2) 

(2-7)' 
27(4 ~ 37) 



(7-2)/2 



(1-7), 



7 



-1/2 



(77) 



and this point is the intersection of the v = surface and 
the magnetosonic critical curve; when h > he, this point 
does not exist. More precisely, one readily finds that for a 
fixed 7, this point moves to infinity as h approaches h~ . 



Expression (77 1 can also be used to determine the loca- 
tion of the 'kink point' of the mEWCSs. In this solution, 
the four forces, viz. thermal gas pressure, magnetic pres- 
sure, magnetic tension and gravity forces are in the same 
order, and the magnetic pressure force is stronger than 
the magnetic tension force with a ratio of —1-1-2/(2 — 7), 
with the magnetic pressure gradient force pointing radi- 
ally outward. In other words, the random magnetic field 
is not force- free. In the isothermal limit of 7 — 1, this 
ratio becomes unity and the random magnetic field is 
essentially quasi-force-free (Yu and Lou, 2005; Yu et al., 
2006). 

Note that this solution can also serve as an asymp- 
totic quasi-static condition, with v approaching zero faster 
than ^ 0{x) and a ~ x^^^". These behaviours are pri- 
marily caused by the polytropic equation of state with 
the magnetic field playing the role of modification. This 
type of asymptotic solutions are referred to as 'quasi- 
static' MHD asymptotic solutions, which can be fur- 
ther sub-divided into two types - type I 'quasi-static' 
MHD asymptotic solutions without oscillations and type 
II 'quasi-static' MHD asymptotic solutions with oscilla- 
tory behaviours. A detailed analysis of this asymptotic 
solution in a polytropic gas without magnetic field has 
been given by Lou and Wang (2006). We here show such 
MHD solutions in our figure illustrations (type II 'quasi- 
static' asymptotic solutions in Figures 17 and 18). An 



analysis of these newly found MHD asymptotic solutions 
is contained in Appendix [G] 
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MHD Solutions Without Crossing the Critical Line 




MHD Solutions Wthout Crossing the Critical Line 



Fig. 7 Global numerical similarity MHD solutions in the 
semi-complete space without crossing the magnetosonic crit- 
ical cur ve c an be rea dily constructed using asymptotic so- 
lutions (39 1 and ( |40[ ) to start at a sufficiently large x (e.g., 
X = 100 j, with the two parameters 7 — 1.01 and n = 0.99 
being specified for a conventional polytropic gas. The two 
perpendicular dotted lines denote the abscissa and ordinate 
axes, respectively, and the straight dotted line —v = —nx is 
the demarcation to the lower left of which solutions become 
unphysical for a negative enclosed gas mass. The two dash- 
dotted curves are the sonic {h — 0) and magnetosonic {h — 1) 
critical curves. The solid {h — 0) and dashed {h = 1) curves 
are similarity MHD flow solutions with /i = and h — 1, re- 
spectively; for the first pair solutions with light linetypes, we 
have both Ao = 4 and Bo = —2; for the second pair solutions 
with heavy linetypes, we have both Ag — 4 and Bg = 0; and 
for the third pair solutions with extremely heavy linetypes, we 
have both Aq — 4 and Bo = 1 as examples of illustrations for 
similarity MHD solutions without crossing the critical curves 
and with outward fiow speeds at far-away region. 
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Fig. 8 The numerical similarity MHD flow solutions with- 
out crossing the magnetosonic critical cur ve c an be read- 
ily constructed using asymptotic solutions ( |39[) and ( |40[ ) to 
start at a sufficiently large x (e.g., x — TDD here), with 
7 = 1.01 and n = 0.99 being specified for a conventional 
polytropic gas. The two perpendicular dotted lines denote 
the abscissa and ordinate axes, respectively, and the straight 
dotted line —v — —nx is the demarcation to the lower left of 
which solutions become unphysical for a negative enclosed gas 
mass. The two dash-dotted lines are the magnetosonic critical 
curves for ft = 10 and h — 100, respectively. The solid and 
dashed curves are MHD solutions with ft — 100 and ft = 10, 
respectively; for the pair solutions with light linetypes, we 
have both Aq — 4 and Bo ~ —15, while for the other pair 
solutions with heavy linetypes, we have both Ao = 4 and 
Bo = —10. The dashed curve with Aq = 4 and Bo = 1 for 
ft = 10 (with an extremely heavy linetype) is an example of 
MHD similarity solutions without crossing the magnetosonic 
critical curve and with an outward fiow speed in the far-away 
region. 



3.4-5 A Global MHD Expansion Solution 



For a constant a in equation (11), we obtain 
v' = 2(x ~ v)/x , 

which can be readily integrated to attain 
u = 2x/3 -|- cx^ , 



(78) 



(79) 



where c is an integration constant. By equation (22 1, we 



should set c = and thus obtain a global MHD solution 
2 



_2x _ 

y ' "~ 3(6/1+1) 



2(n-2/3) , , , 
m= . (80) 



3(6/1-1-1) 

accordingly; the reduced magnetic energy density w is 



w{x) 



9(6/1+1)2 



(81) 



which increases with x quadratically along with the in- 
crease of the radial flow speed v linearly in x. For a poly- 
tropic gas with 7 = 2 — 71, the ratio of the Alfven wave 



speed to the sound speed becomes 

(2-7)/2 





2 




3(6/1 + 1) 



VA 
S 



This solution reaches the singular surface at 

_ r 7(4/3)''/(18/i + 3)^ "1^/^ 

^~ [2{n - 2/3)2/(18/1 + 3) - 4/i/(18/i + 3)^ 

on the magnetosonic critical curve. 



3.4.6 Asymptotic MHD Expansion 

Solutions in the Limit of Large x 



(82) 



Solutions ( 80 ) and ( 81 ) can also be regarded as an asymp- 
totic solution as x approaches infinity. To leading orders, 
the asymptotic MHD solution can be written as 



v{x) = -x + vo , 
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a 



3(6/1+1) 
9(371 + l)(3n 



2)(6/i + 1)2 -6(6/1+1) 



6(n - 3/2)2(3n - 2){<oh + 1) - 2/i(3n - 2) a; 



(83) 



where uq is an integration constant. In this case, the 



denominators of both equations ( 23 1 and ( 24 1 approach 



3(6/1+1) 



2/i 
3(6/i + 



1) 



which does not encounter the niagnetosonic critical curve, 
unless under extremely rare situations. 

3.^.7 Regular Similarity MHD Solutions for Small x 
We may assume an asymptotic MHD series solution as 
a{x) — + aix + a2x'^ + ■ • • , 



v{x) = Uo + ViX + V2x'^ + 



(84) 



as X approaches zero. Substitution of solution (84 1 into 
equations (111 and (22 1 leads to 



v{x) 



157 



(6/i+ l)a* 



.,(2-7) r 



a{x) — a. 



67 



(6/i + l)a* 



w^x) — ha^x^ ~2h- 



67 



(6/i+l)a* 



(85) 



(86) 



(87) 



For a usual polytropic gas with 7 = 2 — n and for very 
small X, the corresponding ratio of the Alfven wave speed 
VA to the sound speed s becomes 



which vanishes as x ^ 0+. As a; approaches zero, the 



denominators of both equations (p3| and (24 1 for this 



series expansion approach —^a^ and these solutions do 
not encounter the magnetosonic singular surface. For this 
asymptotic solution, the four forces including the ther- 
mal pressure, magnetic pressure, magnetic tension and 
gravity forces are in the same order, and the magnetic 
pressure and tension forces tend to be the same in mag- 
nitude. The magnetic pressure gradient force points ra- 
dially outward. 



3.5 Asymptotic Behaviours of Critical Curves 

Asymptotic behaviours of magnetosonic critical curves 
are important from a global perspective. We summarize 
below the major asymptotic behaviours of critical curves 
with or without magnetic field, respectively. 

Case 1 for /i ^ in the presence of magnetic field. 

(i) The limiting regime of a — > +c». According to 
quadratic equation ( 34 1 , the asymptotic quadratic equa- 



tion of X for the magnetosonic critical curve is 



(3n-2)2 



7a 
Jin 



7+1 



2Y 



(89) 



which has only one positive root of small x 

x^ = A{2,n-2fa^-^ , (90) 



and along with equation ( |33[ ), one obtains correspond- 
ingly a diverging radial inflow speed 



(7-i)/2 



(91) 



As a approaches positive infinity, x approaches zero and 
V approaches —00. This appears to be the case also for 
a purely hydrodynamic case with /i = (Lou and Wang, 
2006) and was not discussed in Suto and Silk (1988). 
Thus, the magnetosonic critical curve does not intersect 
the w— axis as compared to the isothermal unmagnetized 
cases (Shu, 1977; Lou and Shen, 2004) where the sonic 
critical line is the straight line —v — \ ^ x. For a magne- 
tized isothermal gas (Yu and Lou, 2005; Yu et al., 2006), 
the magnetosonic critical lines are curves intersecting the 
ti— axis. 

(ii) The limiting case of x — > +00. According to equa- 
tion ( 33 1 , we obtain the following condition for a 



± 



n - 1 + 



a 



{hafl'^ = -n{n-l) 



(92) 



(371-2) 

and correspondingly, the behaviour of v versus x as 

V ^[nT ihay/'^]x . (93) 

This means that as x approaches infinity, v varies linearly 
with X and a approaches a certain constant value deter- 
mined by equation ( |92[ ) . This linear behaviour of u in x is 
qualitatively similar to the isothermal case (Shu, 1977; 
Lou and Shen, 2004), where the sonic critical line is a 
straight line throughout the entire semi-complete space; 
we note that for the isothermal and non-magnetized case, 
a approaches zero as the sonic critical line approaches in- 
finity. For a magnetized gas cloud, the slope of this linear 
relation now depends on h parameter. By the constraint 
of m{x) > 0, the lower signs in solutions (92) and (93 1 



are unphysical, yet mathematically these equations may 
describe asymptotic behaviours of a new branch of crit- 
ical line. 
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MHD Solutions Without Crossing the Critical Line 




Fig. 9 Similarity MHD flow solutions without crossing the 
magnetosonic critical curve with the values of 7 = 1.25 and 
n = 0.75 being specified for a usual polytropic gas. The per- 
pendicular dotted lines denote the abscissa and ordinate axes, 
respectively, and the straight dotted line —v = —nx is a de- 
marcation to the lower left of which the solutions become un- 
physical for a negative enclosed mass. The light dash-dotted 
curves are the magnetosonic critical curves for h = 0, 2, 10. 
The dashed and solid curves are global solutions with h — 
and h = 2, respectively (heavy linetypes: Aq = 4: and Bo = 0; 
light linetypes: Aq = 4 and Bq = —5). The heavy dash-dotted 
curve is the solution with h = 10, Ao = 4 and -Bo ~ —5. The 
dashed curve with Aq — 4 and Bq = 2 for h — (the ex- 
tremely heavy linetype) is an example of hydrodynamic sim- 
ilarity solutions without crossing the sonic critical curve and 
with an outward flow speed in the far-away region. 



The special case of h — he with he defined by equa- 
tion (44 1 should be noted with interest. In this case, the 

-00 becomes 



corresponding value of a as a; 
a = 2(1 - n){3n - 2) , 



(94) 



giving rise to {haY^^ = n and thus the vanishing of the 
leading order term of v by the upper minus sign in equa- 
tion ( 93 1 . This hints that the magnetosonic critical curve 
tends to be parallel to the axis as x approaches infin- 
ity. Naturally for {haY^^ > n as x ^ +00, one infers 
that h > he and vice versa. This reveals an interesting 
trend of variation that as h exceeds the critical value 
he, the magnetosonic critical curve will head up (i.e., 
increase eventually) towards the first quadrant (see sub- 
section 4. 1 ) and h — he marks the marginal case for the 
asymptotic behaviour of the magnetosonic critical curve 
at large x. 

Case II for ft, = in the absence of magnetic field. 
As a — > -l-oo, the asymptotic quadratic equation for 



now becomes 



n^n-lfx^ 



(3n-2)2 



472a2^-2 = , 



(95) 



[see equation (|89]) in parallel for the different coefficient 
of the first term x'*] which has two positive roots for x^. 



(i) For the small x regime of 
= 47(3n - 2)2a'^-3 ^ 



(96) 



one obtains a diverging radial inflow speed of 

As X approaches zero, a approaches positive infinity while 
V approaches negative infinity. Again, the sonic critical 
curve does not intersect the w— axis in contrast to the 
isothermal case (Shu, 1977; Lou and Shen, 2004). Com- 
pared with the isothermal case, this sonic critical curve 
behaviour appears to be a unique feature for a poly- 
tropic gas (Lou and Wang, 2006, 2007), not realized be- 
fore (Suto and Silk, 1988). For an isothermal magnetized 
cloud (Yu and Lou, 2005; Yu et al., 2006), the magne- 
tosonic critical curve involves a cubic equation in terms 
of v; the physical portion of the curve still intersects the 
-y— axis (see figures 1 and 2 of Yu and Lou, 2005). 
(ii) For the large x regime of 



7+1 



n2(n- l)2(3n-2)2 
and accordingly 



(98) 



(99) 



we have the variation trend for the magnetosonic critical 
curve such that as x approaches infinity, a approaches 
infinity and v increases with x linearly. This feature dif- 
fers from the isothermal case (Shu, 1977; Lou and Shen, 
2004); in the isothermal case with n = 1, the above 
asymptotic behaviour is invalid although the asymptotic 
behaviour (99) seems valid. 



3.6 Series Expansions near the 
Magnetosonic Critical Curve 

Generally, smooth analytical solutions cannot go across 



the singular surface as noted in Section 3.3 unless they 
cross the magnetosonic critical curve. Along the mag- 
netosonic critical curve, one cannot calculate the deriva- 
tives of the variables directly from nonlinear MHD ODEs 
( 23 1 and ( 24 1 . From equation (111, we use 



av' — 2a(x — v) jx 
(nx — v) 



(100) 



to compute a' once v' is known. Applying the L'Hopital 
rule to equation (23 1, we immediately obtain the follow- 



ing quadratic equation in terms of v' 

A2(w')'+52W' + C2 = 



(101) 



along the magnetosonic critical curve, where the three 
coefficients A2, B2 and Ci are defined explicitly by 



A. = 1 



(72aT-i -I- 2hax^) 
(ja''^^ + hax"^) 
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MHD Expansion- Wave Collapse Solutions for Different h Values 




MHD Expansion- Wave Collapse Solutions for Different h Values 



Fig. 10 Examples of semi-complete MHD expansion wave 
collapse solutions (mEWCSs) and the corresponding magne- 
tosonic critical curves for different h values with Bq — 0, 
7 = 1.01 and n — 0.99 being specified for a usual polytropic 
gas. The two perpendicular dotted straight lines are the ab- 
scissa and ordinate axes, respectively, and the dash-dotted 
curves are the magnetosonic critical curves. The heavy solid 
curves are the mEWCSs. The corresponding values of Ao for 
the cases of /i = 0, 1, 10 and 30 are 2.0133, 2.0544, 2.5163 
and 5.004, respectively. 



^2 = n + 1 - 4 



(7 a 



2hax'^) 



+4 



C2 



(7a'>'-i 
{-f'^a'-^ + 2hax'^) 



1 



hax^) 
V Ahax 

X 



(nx 



V 



+2 



2(72aT-i + 2hax^ 
(70;'''"^ + hax'^) 



(3n - 2) 
{n - 2)a 



2n 



V 

x"^ 

2hax'^) 



(3n - 2) 



(70^7 1 -|- hax^) 
^{-I'^ai-^ +2hax'^) 



+2ha 



1 



4(x - v) 
{nx — v) 



(102) 



By setting /i = 0, equations (101 1 and ( |102 1 reduce to 
the hydrodynamic results (Suto and Silk, 1988) as a nec- 
essary requirement; by further letting 7 — > 1, these two 
equations reduce to those of the isothermal case (Shu, 
1977). For h ^ and letting 7^1, these two equa- 
tions are equivalent to those of the isothermal MHD 
case (Yu and Lou, 2005; Yu et al., 2006). One can deter- 
mine eigensolution behaviours by Taylor series expan- 
sions near the magnetosonic critical curve using equa- 
tions (101 1 and (1021. By quadratic equation (101 1, one 



obtains two types of eigensolutions across the magne- 
tosonic critical curve. 

For the unmagnetized isothermal case, Whitworth 
and Summers (1985) noticed that there are different be- 
haviours near the vicinity of the sonic critical line, de- 




Fig. 11 Examples of semi-complete mEWCSs and the cor- 
responding magnetosonic critical curves for different h values 
with Bo = 0, 7 = 1.25 and n = 0.75 being specified for a 
usual polytropic gas. The two perpendicular dotted straight 
lines are the abscissa and ordinate axes, respectively, and the 
dash-dotted curves are the magnetosonic critical curves. The 
heavy solid curves represent mEWCSs. The corresponding 
values of Ao parameter for the cases of ^ = 0, 2 and 4 are 
1.151, 2.520 and 21.55, respectively. 



MHD Solutions Without Crossing the Critical Line 



2.520, Bo = 

4, Bo = 
4, Bo = -5 




W 



0.5 



1.5 



2.5 
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X 



Fig. 12 Examples of semi-complete similarity MHD flow so- 
lutions without crossing the magnetosonic critical curve. The 
reduced magnetic energy density w associated with a ran- 
dom transverse magnetic field is displayed in curves (right 
ordinate) using heavy linetypes and the corresponding —v is 
in curves (left ordinate) using light linetypes. The —v curves 
are the same as those in Figure [91 The case of Ao = 2.520 
and Bo = (solid curve) is an mEWCS, while the case of 
Ao = 4 and Bo = (dashed curve; not an mEWCS) is a 
global similarity MHD flow. 
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pending on the topological structure of paths through a 
sonic point. This topological structure is investigated by 
the eigenvalues of the following matrix 



K K 

1 X^- 



(103) 



Table 1 The corresponding m(0) values for similarity MHD 
flow solutions without crossing the magnetosonic critical 
curve for a polytropic gas with 7 = 1.01 and n = 0.99 and 
relatively small values of h parameter. The two coefficients 
Aq and Bo are specified in asymptotic MHD solutions | |39[ ) 
and \40\ . Parameter m(0) is directly computed from equation 
( 10 1. Here, the m(0) value is computed after integrating from 
large x values such as j: = 100 to small x values. 



where A, V and X are the three functions defined by 
equations (251, (26 1 and (27 1, respectively, and Aa de- 



notes a partial differentiation of A with respect to a, 
taking a, v and x as three independent variables; other 
symbols follow the same notational convention by infer- 
ence. The explicit expressions of the partial differentia- 
tions contained in the above matrix are summarized in 
Appendix |B] The characteristic equation for the matrix 
( 103 1 is simply 



A Aa A^ Ax 
-Vo, A-K -Vx 

—Xa —Xy A — Xx 



-0 , 



(104) 



and the signs of the A roots determine the behaviours in 
the vicinity of the magnetosonic critical curve. Equation 
( 104| is equivalent to 

A^ - (A„ + K + Xx)X^ + + Ao^Xx + VyXx 

AxXa AyVa XyVx^X X^AyVx AaVyXx 

+Xo,VyAx + XxVaAy + A^VxXy - V^AxXy = , 

(105) 

which can be cast into a succinct form of 



A^ + B-iX^ + C3A + - , 



(106) 



with apparent definitions for the three coefficients 



C^i and Z?3 by referring to equation (1051. Because the 
magnetosonic critical curve is continuous, this charac- 
teristic equation must have one zero root, i.e., = 0, 
corresponding to a path which stays on the magnetosonic 
curve (see Whitworth and Summers 1985 for a reference 
and Appendix [P] for a pr oof). The other two eigenval- 
ues for A of matrix ( 103 1 are then given by quadratic 
equation 



A^ + S3A + Cg - 



(107) 



The behaviours in the vicinity of the magnetosonic crit- 
ical curve depend upon the signs of the two A roots of 
quadratic equation (107). 

Case I. For a negative determinant B^ — 4C3 < 0, we 
have a spiral or a centre case (e.g., Jordan and Smith, 
1977). In this case, the solutions do not have a one-to-one 
correspondence to x and are thus regarded as unphysical. 
For a polytropic magnetized gas, such points may exist 
(see Appendix |F]) . 

Case II. For a positive determinant B^ — 4C3 > 
with C3 < 0, we have a saddle point along the magne- 
tosonic critical curve with two eigensolutions determined 



h 








1 1 


Ao, Bo 


4, 


4, -2 


4, 4, -2 


m(0) 


5.513 


12.34 


5.912 12.48 



Table 2 The corresponding m(0) values for similarity MHD 
solutions without crossing the magnetosonic critical curve for 
a polytropic gas with 7 = 1.01 and n = 0.99 and relatively 
large values of h parameter. The two coeffi cien ts Ap ari d _Bo 
are specified in asymptotic MHD solutions (391 and (40 1. Pa- 
rameter m(0) is directly computed from equation ( lOrTHere, 
the m(0) value is computed after numerically integrating from 
large x values such as 2; = 100 to small x values. 



h 


10 


10 


100 


100 


Ao, Bo 


4, -10 


4, -15 


4, -10 


4, -15 


m(0) 


40.28 


57.59 


44.58 


59.09 



by equations (100 1, (101 1 and (1021. Numerical integra- 
tions away from these saddle points tend to be stable. 

Case III. For a positive determinant B^ — 4C3 > 
with C3 > 0, we have a nodal point along the mag- 
netosonic critical curve with infinitely many solutions 
crossing the magnetosonic critical curve. As noted by 
Hunter (1986), among these solutions only the two eigen- 
solutions are analytical, while others involve weak dis- 
continuities or weak shocks (e.g., Boily and Lynden-Bell, 
1995) and might be unstable (e.g., Lazarus, 1981). Al- 
though only integrations towards nodal points are stable, 
to pick out the analytic eigensolutions among the so- 
lutions having weak discontinuities, we have integrated 
outward from these points using second-order deriva- 
tives. Explicit expressions of the relevant second-order 
derivatives are summarized in Appendix [C| 

Case IV. For a vanishing determinant i?| — 4C3 — 
0, we have inflection nodal points. If C3 = 0, we have 
degenerate points along the magnetosonic critical curve. 



4 Global Similarity Solutions 

With compatible initial and boundary conditions together 
with a proper treatment of the magnetosonic critical 
curve, the two coupled nonlinear MHD ODEs for self- 
similar collapses and flows can be integrated numerically. 
We have explored MHD solutions numerically, includ- 
ing the properties of the magnetosonic critical curve and 
a — v — X solutions of the M^HD ODEs. In Suto and Silk 
(1988), both cases of n = 1 and n = 2 — 7 were consid- 
ered; here, we focus on the case of a usual polytropic gas 
with n = 2 — 7. In contrast to the case of 7 < 1 consid- 
ered by Fatuzzo et al. (2004), we are mainly concerned 
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Fig. 13 Hydrodynamic polytropic solutions crossing the 
sonic critical curve once, with 7 = 1.01, n = 0.99 and h — 0. 
Point 1 corresponds to a = 5 and x = 0.3941, point 2 corre- 
sponds to a = 3 and x = 0.6570, and point 3 corresponds to 
Q = 2 and x = 0.9877, respectively. The two perpendicular 
light dotted straight lines are abscissa and ordinate axes, re- 
spectively, while the light dotted straight line passing through 
the origin is —v = —nx to the lower left of which solutions 
become unphysical for a negative enclosed mass. The light 
dash-dotted line is the sonic critical curve, and the heavy 
lines are the solutions crossing the sonic critical curve once. 



with the range of 1 < 7 < 4/3. We intend to find semi- 
complete solutions valid in the range of < x < +00, 
and discuss how such MHD solutions can be constructed 
through numerical integrations. 



4.1 Magnetosonic Critical Curves 

The magnetosonic critical curves for different parameters 
can be systematically searched by numerical means, and 
for the completion of a magnetosonic critical curve, one 
needs the relevant analytical results summarized in sub- 



sections 3^ and |3.5| We have extensively explored the 
behaviours of magnetosonic critical curves for specified 
values of n, 7 and h parameters, and present the main 
results below. More details can be found in Appendix [E} 
First, we show magnetosonic critical curves with dif- 
ferent values of h for given n — 0.99 and 7 = 1.01 (Figs. 
[1] and [3]) . A magnetosonic critical curve may be divided 
into two p arts as one picks up different roots of quadratic 
equation ( 34 1 for . 



When h increases for stronger magnetic field strengths, 
the average slope d(— w)/da; of an individual magnetosonic 
critical line increases from negative to positive in our 
figure displays of —v versus x. Meanwhile as the mag- 
netic field becomes strong enough and as x approaches 
infinity, v may approach —00. The critical value he in 



this specific case is 50.521 (see subsection 3.5 1. This fea- 
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Fig. 14 Similarity MHD solutions crossing the magnetosonic 
critical curve once, with 7 = 1.01, n — 0.99 and /i = 1 be- 
ing specified for a conventional polytropic gas. Point 1 cor- 
responds to a = 15 and x = 0.1200, point 2 corresponds 
to a = 5 and x — 0.3214, point 3 corresponds to a = 2 and 
X = 0.7021, and point 4 corresponds to a = 1 and x = 1.2456, 
respectively. The two perpendicular light dotted straight lines 
are abscissa and ordinate axes, respectively, while the light 
dotted straight line passing through the origin is —v = —nx 
to the lower left of which solutions become unphysical for a 
negative enclosed mass. The light dash-dotted line is the mag- 
netosonic critical curve, and the heavy lines are the solutions 
crossing the magnetosonic critical curve once. 



MHD flow solutions not crossing the magnetosonic crit- 
ical curve. In contrast to the straight sonic critical line 
for the isothermal unmagnetized case (Shu, 1977; Lou 
and Shen, 2004), the magnetosonic critical curves here 
diverge as x approaches zero (see Fig. |3]). 

The magnetosonic critical curves in this nearly isother- 
mal case of 7 > 1 can be compared with those of the 
isothermal MHD case of 7 = 1 (Yu and Lou, 2005; Yu 
et al., 2006). Their asymptotic behaviours are different 
for both limiting regimes of x — *■ 0^ and x — > +00. 
As X approaches zero, the magnetosonic critical curve in 
the isothermal case intersects with the vertical t)— axis, 
while in the nearly isothermal MHD case it diverges as 
X ^ 0^ . As X approaches infinity, the magnetosonic crit- 
ical curve in the isothermal case remains in the fourth 
quadrant, while in the nearly isothermal case it can head 
up to the first quadrant. These quaUtative differences in 
asymptotic behaviours in such parallel cases result from 
equations (91 1 — (93). According to equation (91 1, for 
V remains finite as x approaches zero, while for 



ture is important in the numerical analysis of similarity 



7=1, 

7 > 1, even a small increment in 7 will lead to a diver- 
gence of V as X goes to zero. In accordance with equa- 
tion (92 1, if n = 1 as in the isothermal case, a must 
approach zero as x approaching infinity, which means 
that V nx, i.e. the magnetosonic critical curve remains 
in the fourth quadrant. Nonetheless when n is not equal 
to unity, this constraint on asymptotic behaviour of the 
magnetosonic critical curve disappears. Another perspec- 
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Table 3 The corresponding m(0) values for self-similar 
MHD flow solutions without crossing the magnetosonic crit- 
ical curve for a conventional polytropic gas of 7 = 1.25 and 
n = 0.75 with different values of h. The two param eter s Aq 
and Bo are specified in asymptotic MHD solutions ( 39 1 and 
(40 1. Parameter m(0) is directly computed from equation 
(TDl. Here, the m(0) value is computed after integrating from 
large x values such as a; = 100 to small x values. 



MHD Solutions Crossing the Critical Line Once 



h 








2 


2 10 


Ao, Bo 


4, 


4, -5 


4, 


4, -5 4, -5 


m(0) 


3.8185 


4.952 


3.664 


4.919 4.747 



Table 4 For Figure 10 we summarize the corresponding 
m(0) values and the x values of the 'kink point' x^ for MHD 
expansion wave collapse solutions for 7 — 1.01 and n — 0.99 
for a conventional polytropic gas (see Yu & Lou 2005 and Yu 
et al. 2006 for an isothermal magnetized gas). 



h 





1 


10 


30 


m(0) 


1.0120 
1.019 


1.4931 
1.76 


4.346 
5.08 


19.26 
12.11 




Table 5 For Figure |ll[ we summarize the corresponding 
m(0) and the x values of the 'kink point' values Xk for 
MHD expansion wave collapse solutions with two parame- 
ters 7 = 1.25 and n — 0.75 specified for a usual polytropic 
gas (see Yu and Lou 2005 for an isothermal magnetized gas). 



h 





2 


4 


m(0) 


0.7992 


2.083 


24.38 


Xk 


1.37 


2.52 


6.68 



five is that when n — s- 1, we have /ic 00; thus whatever 
h values will lead to asymptotic behaviours such that 
V — > +00 as X +00. 

The magnetosonic critical curves for different values 
of h given n = 0.75 and 7 = 1.25 are shown in Fig. 
[2] When h increases, again the average slope d{—v)/dx 
of an individual magnetosonic critical line increases from 
negative to positive in the —v versus x presentation. The 
value of he is 4.5 in this example. As the magnetic field 
becomes extremely strong, one obtains another branch 
of the magnetosonic critical curve as shown in Fig. |2]for 
h = 100. This new branch is the one mentioned in equa- 
tion ( 33 1 . The lower branch of the heavy dotted curve is 



unphysical for being to the lower left of the straight line 
—V = —nx. Also the magnetosonic critical curve diverges 
as X approaches zero. 

Enlarged features for diverging behaviours of v{x) 
along the magnetosonic critical curves for small x in Figs. 
[l]and|2]are shown in Figs. [3] and |4] respectively. The cor- 
responding features of a versus x are displayed in Figs. [5] 
and [6] respectively. The basic facts that for ft, = 0, a ap- 
proaches infinity both as x approaches zero and infinity, 
while for h > 0, a approaches infinity as x approaches 
zero and a approaches a constant as x goes to infinity 
are all consistent with the relevant analytical results pre- 
sented in subsection 13.51 



(a) The smaller x regime. Crossing points 1, 2, and 3 corre- 
spond to a = 5 and x — 0.2170, a — 0.9 and x — 0.8255, and 
a = 0.3 and x — 1.9389, respectively. 

MHD Solutions Crossing the Critical Line Once 



> 

I 




X 

(b) The larger x regime. Crossing points 1, 2, and 3 corre- 
spond to a = 0.1 and x = 4.8147, a = 0.05 and x = 9.9587, 
and a — 0.04 and x — 14.0005, respectively. 

Fig. 15 MHD solutions crossing the magnetosonic critical 
curve once, with 7 = 1.01, n = 0.99 and h = 10. Panel 
(a) displays the solutions crossing the magnetosonic critical 
curve at smaller x, while panel (b) displays those crossing the 
magnetosonic critical curve at larger x. The two perpendicu- 
lar light dotted lines are abscissa and ordinate axes, respec- 
tively, and the light dotted straight line is for —v = —nx. 
The light dash-dotted line is the magnetosonic critical curve, 
and the heavy curves are the MHD solutions crossing the 
magnetosonic critical curve once. 



By numerical exploration, we found that the magne- 
tosonic critical curve has two branches in the 7 = 1.01, 
n — 0.99 case with h = 1000. Also, the critical curve con- 
sists of two parts as shown in Fig. [ijin the case oi h = 
can also be found in the case of 7 = 1.01, n = 0.99, and 
h = 0.03. From the results for critical curves, one can 
see that asymptotic analyses in subsections |3.3| and |3.5| 
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Fig. 16 MHD solutions crossing the magnetosonic critical 
curve once with parameters 7 = 1.01, n = 0.99 and h = 100 
being specified. Crossing point 1 corresponds to a = 0.2 and 
X = 1.5055, point 2 corresponds to a = 0.05 and x = 4.8487, 
and point 3 corresponds to a = 0.024 and x = 12.2494, re- 
spectively. The two perpendicular light dotted straight lines 
are abscissa and ordinate axes, respectively, and the light 
dotted straight line is the demarcation line —v = —nx. The 
light dash-dotted line is the magnetosonic critical curve, and 
the heavy curves are the solutions crossing the magnetosonic 
critical curve once. 



are necessary for determining the entire magnetosonic 
critical curve. 



4.2 Solutions without crossing the Magnetosonic 
Critical Curve and MHD Expansion Wave Collapse 
Solutions 

4-2.1 General MHD Solutions without Crossing the 
Magnetosonic Critical Curve 

Among the asymptotic MHD solutions derived in sub- 
section [3]4|the series expansion at large x described by 
equations (|39| and (40) can be readily integrated from 



large values of x inward to obtain numerical solutions 
without encountering the magnetosonic critical curve. 
Specifically, we integrate the solutions from a starting 
point oi X = 100. 

We present such global semi-complete similarity MHD 
solutions for the case of 7 = 1.01 and n — 0.99 in both 
Figures |7] and |8] Note that when x approaches zero, the 
solutio n app roaches the free-fall state as discussed in sub- 
section |3!4] [see equations (50 1 and (51)]. Note also that 
the major difference in v{x) of the two solutions with 
the same values of and -Bo but with different mag- 
netic field strengths (i.e., different h) manifests mainly 
at small x about < x < 10 in both Figures. Both Figs. 
[7] and [8] show that the magnetic field mainly accelerates 
the central collapses, i.e. when h is larger, the —v{x) 
becomes larger at the same x, although the asymptotic 
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Fig. 17 Hydrodynamic solutions crossing the sonic critical 
curve once, with 7 = 1.25, n — 0.75 and h = 0. Crossing 
point 1 corresponds to a — 5, x = 0.1386, crossing point 2 
corresponds to a = 2 and x = 0.3166, and crossing point 
3 corresponds to a = 1 and x = 0.6131, respectively. The 
two perpendicular light dotted lines are the abscissa and or- 
dinate axes, respectively, and the light dotted straight line 
—V = —nx is the demarcation line to the lower left of which 
solutions become unphysical. The light dash-dotted line is 
the sonic critical curve, and the heavy curves are the solu- 
tions crossing the sonic critical curve once. The oscillatory 
behaviours of solutions Hi, 2Ii and 3Ii as x approaches are 
not readily seen and will be discussed at the end of the main 
text. 



behaviours of v(x) as x approaches infinity show that 
larger h implies smaller —v{x) at the same large x. 

Numerical similarity MHD solutions for the case of 
7 — 1.25 and n — 0.75 are presented in Figure [9] In this 
figure the solution with the same values of Aq and Bq but 
larger h value cannot catch up with the other solutions 
with smaller values of h. This does not necessarily mean 
that the net magnetic force does not accelerate collapses, 
because of a smaller in the initial state for the case 
of a larger h. 

We briefly note several points here. Firstly, we ana- 



lyzed in subsection 3.4 the asymptotic behaviour of the 
MHD free-fall solutions for small x and inferred from 
that analysis that these solutions will not encounter the 
magnetosonic singular surface at small x, meanwhile the 
asymptotic MHD solutions as x approaches infinity [equa- 
tions (39) and (40)] also do not encounter the singu- 
lar surface. Therefore, these numerical MHD solutions 
are specific examples of semi-complete solutions without 
crossing the magnetosonic critical curve or encounter- 
ing the singular surface. Secondly, there exists a two- 
dimensional continuum regime of parameters Aq and Bq 
for this series expansion of MHD solution, e.g., for a spec- 
ified Aq parameter, parameter Bq should be larger than a 
certain threshold value in order to construct MHD sim- 
ilarity fiow solutions without encountering the magne- 
tosonic critical curve. Outside such allowed parameter 
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Fig. 18 MHD solutions crossing the magnetosonic critical 
curve once with 7 = 1.25, n = 0.75 and h — 2 being speci- 
fied. Crossing point 1 corresponds to a = 1 and x = 0.5107, 
crossing point 2 corresponds to a = 0.4 and x — 1.1549, and 
crossing point 3 corresponds to q = 0.18 and x — 3.7395, 
respectively. The two perpendicular light dotted lines are the 
abscissa and ordinate axes, respectively, and the light dot- 
ted straight line —v = —nx is the demarcation to the lower 
left of which solutions become unphysical. The light dash- 
dotted curve is the magnetosonic critical curve, and the heavy 
curves are the solutions crossing the magnetosonic critical 
curve once. The oscillatory behaviours of solutions Hi, 2Ii 
and 3Ii as x approaches are not readily seen and will be 
discussed at the end of the main text. 



regime, the solutions tend to crash on to the singular 
surface but away from the magnetosonic critical curve 
so that a global semi-complete MHD solution does not 
exist. Thirdly, the MHD similarity flow solutions with 
h > he will cross the magnetosonic critical curve at the 
projection to —v ^ x plane, yet they do not actually 
cross the magnetosonic critical curve in the a ~ v ~ x 
space because they do not encounter the singular sur- 
face. 



4-2.2 Construction of MHD Expansion 

Wave Collapse Solutions (mEWCSs) 

One interesting solution among global MHD solutions 
not crossing the magnetosonic critical curve is the lim- 
iting solution corresponding to MHD expansion wave 
collapse solutions (mEWCSs) as a generalization of the 
isothermal EWCS (Shu, 1977) in two important aspects, 
i.e., the poly tropic gas and the inclusion of a random 
magnetic field. We emphasize here the existence of such 
global MHD similarity solutions because magnetic fields 
do exist in molecular clouds in general and play impor- 
tant roles in the evolution of a collapsing cloud. From the 
perspectives of dynamic evolution, diffusive processes, 
radiative signatures, formations of discs and jets and 
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Fig. 19 MHD solutions crossing the magnetosonic critical 
curve once with 7 = 1.25, n = 0.75 and h = 10. Crossing 
point 1 corresponds to a = 1 and x — 0.3998, crossing point 
2 corresponds to a = 0.26 and x — 1.2053, and crossing 
point 3 corresponds to a = 0.12 and x = 4.6518, respectively. 
The two perpendicular light dotted lines are the abscissa and 
ordinate axes, respectively. The straight line —v = —nx is a 
demarcation line to the lower left of which solutions become 
unphysical. The light dash-dotted curve is the magnetosonic 
critical curve, and the heavy curves are the MHD solutions 
crossing the magnetosonic critical curve once. 
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Fig. 20 Correspondin g a be haviours for the MHD solutions 
IH, 211 and 3H in Fig. 15(a) The light dash-dotted curve is 
the magnetosonic critical curve, and the heavy curves are the 
solutions crossing the magnetosonic critical curve once. 



origin of stellar magnetic fields, one must take into ac- 
count of magnetic fields in molecular clouds. By the 
model scenario here, our analysis suggests that there ex- 
ist mEWCSs for h < he {a. weaker magnetic field), while 
no such mEWCS exists for h > he (a stronger magnetic 
field). 

This conclusion can be viewed from the following 



perspectives. First, in asymptotic MHD solutions (39 1 
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and ( 40 1 , parameter Bq can be regarded as an exter- 



Phase Diagram for Hunter Type Solutions 



nal (or initial) radial flow speed more or less indepen- 
dent of the mass density profile. For example in the 
isothermal case, parameter Bq represents an asymptotic 
steady flow speed in regions fa r fro m the core (see Lou 
and Shen 2004 and subsection 3.4 here). Parameter Aq 
contributes to the radial speed profile due to gas mass 
density distribution and the associated self-gravity. To 
construct mEWCSs, one should require that Bq = and 
the reduced radial speed v{x) approaches 0~, i.e., for the 
limiting series solution, the radial fiow velocity remains 
negative and approaches zero far away. In this limiting 
regime, a magnetized gas cloud of quasi-spherical sym- 
metry remains at rest in early times and the core collapse 
is induced by self-gravity. In this perspective, the critical 
Aq for mEWCS is determined by setting the coeflicient 
of x'^""^-'/" term [from the same two terms x'""^'/" and 
2;(2-27-n)/ri ^ usual polytropic gas with n — 2 ~ j] 
in asymptotic MHD solution (40 1 to vanish. In contrast, 
the case oi h > he corresponds to positive coeflicients 
of both terms involving Aq in asymptotic MHD solution 



(40) because of l/(3n -2) + 2h{n - l)/n^ < 0, indicat- 



ing outward expansions for any value of Aq > when 
the coeflicient of x*^""^^/" term vanishes. Secondly, for a 
mEWCS, the reduced velocity remains zero at large x 
until the solution meets the magnetosonic critical curve 
at the point where the magnetosonic critical curve inter- 
sects the cc— axis. At this intersection point, the slope of 
the magnetosonic critical curve is negative for h < he, 
which allows the mEWCS to head up as a; — > O"*", while 
for h > he, the slope of the magnetosonic critical curve 
is positive to force the solution to crash on to the sin- 
gular surface without leading to a mEWCS. Thirdly, the 
physical reason for the non-existence of a mEWCS when 
h > hcis that a sufficiently strong magnetic field tends to 
prevent a gravitational collapse and to drive an outward 
expansion instead. The stability analysis of the present 
MHD problem remains to be examined for the case of 
h > he- For h > he a, magnetized gas cloud collapses 
only when external (initial) inflows are present. Because 
of he oo in the isothermal case, any h may lead to 
an mEWCS. Finally, when h > he or Aq smaller than a 
certain critical value representing the mEWCS, parame- 
ter Bq should be negative to ensure the core collapse of 
a magnetized gas cloud. When Bq becomes sufficiently 
negative, the corresponding initial flow speed proflle will 
have a tendency to collapse, and a core collapse without 
encountering the magnetosonic critical curve can hap- 
pen. Therefore for any Aq > 0, one can find solutions 
without encountering the magnetosonic critical curve. 



In Figures nO\ and 11 



we present the mEWCSs of 
a usual polytropic gas cloud for 7 = 1.01, n — 0.99 
and 7 = 1.25, n = 0.75, respectively. We conclude from 
both Figures 10 and [TT] that the kink points Xk of the 
mEWCSs have larger x values as h increases. The corre- 
sponding values of m(0) and Xk are summarized in Tables 
|4]and|5] respectively. 
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Fig. 21 Density-speed phase diagram for MHD generaliza- 
tions of Hunter type solutions with 7 = 1.01, n = 0.99, h — 1, 
and a chosen meeting point at xp — 0.5. The 'outward' curve 
represents the phase path of Hunter type solutions when the 
parameter a, is changed gradually. The 'inward' curve repre- 
sents the phase path of solutions crossing the magnetosonic 
critical curve as the intersection point of the MHD solutions 
with the magnetosonic critical curve is gradually adjusted. 



We also present the reduced magnetic energy density 
w = ha^x^ associated with a random transverse mag- 
netic field Bf versus x. Figure [12] collects a sample of 
w versus x solution curves for 7 — 1.25, n — 0.75 and 
h = 2 (see Fig. [9|. Note that the initial magnetic field 
strengths are the same as x approaches in- 

finity (i.e., t — !■ 0"*") for the two upper curves (i.e., dashed 
and dash-dotted linetypes) with Aq = 4; at first the 
magnetic field strength increases faster in the case of 
Bq = 'AS X becomes smaller, and then at some point 
around x ~ 0.63549 the curve with larger initial veloc- 
ity {Bq ~ —5) catches up with the former and grows 
faster. For the mEWCS with Aq = 2.520 and Bq ^ 
and a decreasing x, the reduced magnetic energy den- 
sity w begins to increase before the kink point Xk = 2.6 
in the reduced velocity field, and at the kink point the 
magnetic field also appears to increase more slowly than 
at larger x values. The interpretation of this feature is 
that at a specific point, the magnetic field first main- 
tains a constant and then begins to decrease when the 
magnetosonic wave front reaches this point. 

One readily obtains the corresponding value of m(0) 
for each MHD free-fall solution using equation ( [To| ) . The 
corresponding m(0) values for all solutions computed in 
this section are summarized in Tables [T| to [5j According 
to Tables [l] and |2] one finds that m(0) increases with 
either larger h or larger magnitude of Bq, indicating that 
a stronger magnetic field and a faster inward initial speed 
both result in a more rapid core collapse and lead to an 
enhanced central mass accretion (n.b. profiles of gas mass 
density scalings remain the same) . This is related to the 
fact that the inward magnetic tension force is twice as 
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Fig. 22 The first three discrete MHD Hunter type solutions 
for 7 = 1.01, n = 0.99, /i = 1 and xp = 0.5 in the semi- 
complete space. The two perpendicular light dotted lines are 
the abscissa and ordinate axes. The light dash-dotted curve 
is the magnetosonic critical curve and the heavy lines are 
the similarity MHD solutions curves. Curve 1 corresponds to 
Q, — 359.37 and crosses the magnetosonic critical curve at 
X = 2.1154, curve 2 corresponds to a* = 8.3941 x lO'' and 
crosses the magnetosonic critical curve at a; = 1.7248, and 
curve 3 corresponds to a* — 1.4769 x 10^^ and crosses the 
magnetosonic critical curve at a; = 1.7637, respectively. As 
X — > -l-oo, the two corresponding Ao and So in asymptotic 
solutions (l39| and (|40| are Aq = 2.430 and Bo = 0.3810 for 
curve 1, 2.0lVand Bo = -0.03998 for curve 2, and 

Ao = 2.058 and Bq — 0.003869 for curve 3, respectively. 

Tabl e 6 Values of the two parameters Aq and Bo in equa- 
tions ( 39 1 and ( 40 1 at large x of relev ant semi-co mpl ete sim- 
ilarity MHD flow solutions in Figures [TS] through |19[ 
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Table 7 Corresponding m(0) values of MHD free-fall solu- 
tions at small x for relevant semi-complete solutions in Figs. 
[T3l through PI 
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Fig. 23 Enlarged versions for the MHD generalizations of 
the first three Hunter type MHD solutions for 7 = 1.01, n = 
0.99, h = 1 and xp = 0.5 to illustrate that the ith solution has 
i stagnation points (i = 1, 2, 3). Other parameters are the 
same as those in Figure [22] Note that the ordinate scale has 
a factor of 10~* and the abscissa for x is in the logarithmic 
scale. The feature of self-similar magnetosonic oscillations is 
shown by these MHD solutions, as a general extension of the 
isothermal hydrodynamic feature revealed by Hunter 1977 
and Lou & Shen 2004. 



large as the outward magnetic pressure gradient force. In 
reference to Table |4] we note the increase of m(0) with 
a larger h in niEWCSs with the two parameters 7 and 
n being specified. By Table [3j one finds that a larger 
Bq results in a smaller m(0), but as mentioned above 
this does not necessarily mean that the magnetic field 
does not accelerate the core collapse because of an initial 
smaller velocity. By Table |5] one again gets an increase 
of m(0) with increasing h. We also provide the x values 
of the corresponding 'kink points' of the mEWCSs 
in Tables |4] and [5] showing that for a given 7 > 1 , the 
corresponding Xk increases with an increasing h value as 
long as h < he- 



4.3 Numerical MHD Solutions Crossing the 
Magnetosonic Critical Curve Once 

Behaviours of similarity MHD fiow solutions around the 
magnetosonic critical curve are determined by equations 
(100 1 — (102), where the first derivatives of v{x) and 



a{x} with respect to x can be determined along the mag- 
netosonic critical curve. Numerically, one can integrate 
from a point in the vicinity of the magnetosonic critical 
curve away to obtain a portion of the eigensolution cross- 
ing the magnetosonic critical curve. At one specific point 
on the magnetosonic critical curve, there exist two eigen- 
solutions crossing the magnetosonic critical line. The one 
of smaller v' in the vicinity of the magnetosonic criti- 
cal line is denoted as type I solutions and the other of 
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Phase Diagram for Xp = 0.5 at 7 = 1.01 and h = 




outward 



-0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0.1 0.2 0.3 

V 

Fig. 24 Density-speed phase diagram of a versus v for 7 = 
1.01, n = 0.99, h — and xp = 0.5. The 'inward curve' 
(dashed curve) denotes the phase path by integrating inward 
from different x {x > xp) along the sonic critical curve to 
reach xp, using the type I eigensolution, while the 'outward 
curve' denotes the phase path integrating from different x 
[x < Xp) along the sonic critical curve to reach xp, using 
the type II eigensolution. The spiral pattern of the outward 
curve suggests the trend of infinitely many matches, leading 
to infinitely many semi-complete polytropic EECC solutions 
(see Lou and Shen 2004). 



larger v' is type II solutions [for a sufficiently small x, 
these parallel with type 1 and type 2 derivatives, respec- 
tively (Lou and Shen, 2004)]. These notations differ from 
those used in the isothermal case (Shu, 1977; Lou and 
Shen, 2004). We have explored MHD solutions crossing 
the magnetosonic critical curve once, with several typi- 
cal results shown in Figs. [13] to |19[ For mnemonics, we 
denote the type Y (Y=I, II) solution outward (inward) 
from the xth {x — 1, 2, 3, 4) point by a;Yo (xYi) solu- 
tion. We searched similarity MHD solutions for the cases 
of 7 = 1.01, n = 0.99 and of 7 = 1.25, n = 0.75 with 
various h values, i.e., different reduced magnetic energy 
density. We mainly focus on semi-complete physical so- 
lutions in 0+ < a; < +00. 



4-3.1 Solutions with 7 = 1.01 and Small h Values 

Hydrodynamic solutions with 7 = 1.01, n — 0.99 and 
h = are shown in Figure 13 Solutions lli, 2Ii and 3Ii 
for small x values all run under the straight demarca- 
tion line — w = —nx (to the lower left of which solutions 
become unphysical) and then encounter the singular sur- 
face (without crossing the sonic critical curve); they are 
thus not valid at small x neither mathematically nor 
physically. Solutions lllo, 2IIo and 3IIo crash on to the 
singular surface and are invalid at large x mathemati- 
cally, although one may expect special solutions crossing 
the critical line twice at some specific points in a discrete 
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Fig. 25 Corresponding to the matches in the density-speed 
phase diagram of Figure |24[ the first four hydrodynamic 
polytropic solutions with 7 — 1.01, n — 0.99, h — and 
Xp — 0.5 are shown here. Curve 1 (solid line) is obtained 
from X = 1.7248 {a = 1.15739, v = 0.7018) and x = 0.2035 
{a = 9.7036, v = —0.8150) along the sonic critical curve 
(this is a type I— type II match), curve 2 (dashed line) is ob- 
tained from X = 0.7448 (a = 2.64742, v = -0.2726) and 
X = 2.2124 X lO"" (a = 9224.0, v = -1.0517) along the 
sonic critical curve, curve 3 (dotted line) is obtained from 
X = 1.1319 (a = 1.74780, v = 0.1128) and x = 4.9114 x lO"*^ 
{a = 4.2354xl0^ v = -1.0722) along the sonic critical curve, 
and curve 4 (dash-dotted line) is obtained from x = 0.9855 
(a = 2.00440, v = -0.0329) and x = 2.8282 x 10~* (a = 
7.5482 X 10^ V = -1.1004) along the sonic critical curve. 
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Fig. 26 Corresponding to Figs. 24 and 25 the enlarged ver- 
sions of the first four solutions for smaller a:: in a logarithmic 
scale with 7 = 1.01, n = 0.99, h — and xp = 0.5 are shown 
with stagnation points and radial oscillations. Curves 1, 2, 3 
and 4 are the same as those in Fig. |25| 



manner (see subsection 4.4 and Lou and Shen 2004. So- 
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Phase Diagram for Xp = 0.5 at 7 = 1.01 and h = I 
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Fig. 27 Density-speed phase diagram for specified parame- 
ters 7 = 1.01, n = 0.99, h = \ and a chosen meeting point 
Xp = 0.5. The 'inward curve' is the phase path obtained by 
integrating inward from different x {x > xp) along the mag- 
netosonic critical curve to reach xp, using the type I MHD 
eigensolution, while the 'outward curve' denotes the phase 
path integrating from different x {x < xp) along the mag- 
netosonic critical curve to reach xp, using the type 11 MHD 
eigensolution. The spiral pattern of the 'outward curve' hints 
at infinitely many matches, leading to infinitely many semi- 
complete MHD EECC similarity solutions. 
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Fig. 28 Corresponding to the three meeting point matches 
in the phase diagram of Figure |27| we show the first three 
MHD similarity solutions with 7 = 1.01, n = 0.99, h — 1 and 
Xp — 0.5. Curve 1 (solid line) is obtained from x — 2.1065 
(a = 0.526485, v = 0.2580) and x = 0.001157 (q = 1747.8, 

V = —1.0432) along the magnetosonic critical curve, curve 
2 (dashed line) is obtained from x = 1.7248 {a = 0.671990, 

V = -0.02598) and x = 5.5880 x 10"** (a = 3.8072 x 10^ 

V = —1.0966) along the magnetosonic critical curve, and 
curve 3 (dash-dotted line) is obtained from x = 1.7637 
(a = 0.653929, v = 0.002533) and x = 3.3391 x 10"^^ 
[a = 6.6905 x lO^'^, v = —1.1516) along the magnetosonic 
critical curve, respectively. For the convenience of visual in- 
spection, we have multiplied curve 3 by a factor of 5, i.e., ~5v 
is shown here for curve 3. 
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Fig. 29 The enlarged version of the first three MHD EECC 
solutions with 7 = 1.01, n = 0.99, h = 1 and xp = 0.5 in a 
logarithmic scale for small x to show the stagnation points 
with V — and self-similar magnetosonic radial o scill ations. 
Curves 1, 2 and 3 are the same as those in Figure [28] 



lutions llli, 2IIi and 3IIi approach the free-fall solution 
[equations (50 1 and (51 1] and are valid for small x. Solu- 
tions llo, 2I0 and 3Io approach asymptotic solutions ( 39 1 



and (40 1 and are also valid at large x. Solutions shown in 
Fig. |13| fail to form semi-complete hydrodynamic poly- 
tropic solutions. Likewise, MHD solutions with the same 
7 and n but h — 1 shown in Fig. [14] also fail to form 
semi-complete MHD solutions, and the validity of each 
MHD solution is just like the corresponding one in Fig. 
[T3I i.e., MHD solutions Hi, 2Ii and 3Ii run to the lower 
left of the demarcation line — w = — nx and then crash 
on to the singular surface, solutions lllo, 2IIo and 3no 
crash on to the singular surface, solutions IHi, 2ni and 
3ni approach the MHD free-fall solution and solutions 
llo, 2I0 and 3Io approach asymptotic MHD solutions 
dsgl and (|40|. 



4-3.2 Solutions with 7 = 1.01 and Larger h Values 

For a larger h value of stronger magnetic field, we show 
a few more examples of semi-complete MHD flow solu- 
tions. Fig.[T5|shows the solutions with 7 — 1.01, n — 0.99 
and h = 10. T he solutions lli, 2Ii, 3Ii in Fig. 15(a) and lli 
in Fi g. |15(b) approach asymptotic MHD solutions (62 1 
" ;h C = 7.62 which is the same for all four 



and ([63| wit^ 

solutions here. Solutions Ho, 2Io and 3Io in both pan- 
els (a) and (b) approach the asymptotic forms of (39 1 
and (40 1 at large x, with the corresponding parameters 
Aq and Bq summarized in Table [6j Solutions llli, 2IIi 
and 3IIi in both panels (a) and (b) approach the MHD 
free-fall solution at smaller x with to(0) summarized in 
Table It] Solutions lllo, 2IIo and 3IIo in both panels (a) 
and (b) approach asymptotic MHD solution ( 83 ) . Solu- 
tions 2Ii and 3Ii in Fig. ]15(b)] run to the lower left of the 
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Phase Diagj'am for xp = 0.2 at 7 = 1.25 and h 
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Fig. 30 Density-speed piiase diagram for 7 = 1.25, n = 0.75, 
h = and xp = 0.2. The 'inward curve' is tiie phase curve 
obtained by integrating inward from different x {x > xp) 
along the sonic critical curve to reach the chosen xp using 
the type I eigensolution, and the 'outward curve' denotes the 
phase path integrating from different x (x < xf) along the 
sonic critical curve to reach xp using the type II eigensolu- 
tion. There is no spiral pattern for the 'outward curve' and 
only one single match can be found at the lower-right corner. 
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MHD Solutions Crossing the Critical Line Twice 
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Fig. 31 Corresponding to the single match in the density- 
speed phase diagram of Figure [30l the relevant hydrodynamic 
polytropic solution —v{x) versus x (solid curve) is shown here 
for 7 = 1.25, n = 0.75, h = Q and xp = 0.2, with the two 
crossing points given hy x = 3.1529 (a = 0.36876104, v — 
1.3777) and x = 0.01197 (a = 80.978, v = -1.9275) along 
the sonic critical curve (dash-dotted curve). 



demarcation line —v ~ —nx and then crash on to the 
singular surface (not crossing the magnetosonic critical 
curve) and are thus invalid in the small x regime. There- 
fore in Fig. [15] among the twelve solution branches (II, 
12, 13, III, 112 and 113 in both panels), there exist in to- 
tal ten semi-complete solutions, all but the two solution 
branches 21 and 31 in Fig. |15(b)] In Fig. 16 the solutions 



are also all semi-complete solutions and the correspond- 
ing parameters are tabulated in Tables [6| an d [T} The cor- 
responding C parameter in equations ( |62[ ) and ( |63[ ) for 
this case is C = 95.0. The novel asymptotic MHD solu- 



tions ( 62 1 and ( 63 1 have been used here and the relevant 



numerical MHD results are consistent with the analytical 
analysis. 

4-.3.3 Solutions with 7 — 1.25 and Small h Values 

For a usual polytropic gas with a different combination 
of larger 7 and smaller n, e.g., 7 = 1.25 and n = 0.75, the 
situation remains qualitatively similar, but in the smaller 
X regime with also a small h, the situation appears some- 
what different. Fig. |17| presents the unmagnetized case 
of ft, = where solutions lli, 2Ii and 3Ii show self-similar 
oscillation behaviours about the x axis for the regime of 
small x (see Fig. 36 and the isothermal case Lou & Shen 
2004). This MHD oscillation behaviour, however, differs 
from that described by Lou and Shen (2004) , mainly be- 
cause it is of type II 'quasi static' (i.e., the infall speed 
remains finite) asymptotic solution as was recently dis- 
covered in the hydrodynamic case (Lou and Wang, 2006); 
the presence of such 'quasi-static' asymptotic solutions 
is intimately related to the polytropic equation of state. 
Given other parameters the same but for h — 2, Fig. [18] 
displays semi-complete MHD solutions III, 211 and 311, 
with the inner portions approaching the MHD free-fall 
state with parameter m(0) summarized in Table |7] and 
the outer portions approaching asymptotic MHD solu- 
tions (83 1. Meanwhile, MHD solutions lli, 2Ii and 3Ii 



clearly display the self-similar magnetosonic oscillation 
behaviours for smaller x in an analogous manner. This 
confirms the existence of the 'quasi-static' asymptotic 
behaviour in our MHD model. An analysis of such MHD 
solutions is contained in Appendix [G} 

4-.3.4 MHD Solutions with 7 ~ 1.25 and Larger h Values 

For an even larger h such a,s h ~ 10, one can construct 
semi-complete MHD solutions for both types. Besides the 
similar type of III, 211 and 311 solutions with parameter 
to(0) also contained in Table [t] MHD solutions II, 21 
and 31 are also valid both in small and large x regimes. 
The corresponding C in solutions ( [62| ) and ( 63 1 for MHD 
solutions Hi, 2Ii and 3Ii is C = 1.47, and the parameters 
Aq and Bq of equations dMl) and (|40l) are summarized 



in Table [6j The self-similar magnetosonic oscillation be- 
haviour is described in section |5] The novel asymptotic 
MHD solutions also exist for 7 = 1.25. 

4-3.5 Solution Behaviours of a versus x 

Behaviours of the corresponding a{x) so lution profile of 
the MHD EECC solutions in Fig. 15(a) are presented in 
Fig. [20] Note that the a{x) profiles of all three MHD solu- 
tions approach the same limiting constant 2/[3(6/i-M)] at 
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large x as expected. Other MHD EECC solutions in the 
figures above are similar to these solutions. This type of 
solutions represents a magnetized gas in spherical enve- 
lope expansion with a constant reduced density a at large 
X (initially) and a free-falling core for small x (finally). 
For a specified set of 7, n and h, this series of MHD 
solutions may form a one-dimensional continuum, de- 
pending on the point where they cross the magnetosonic 
critical curve. In other words, by adjusting the cross- 
ing point where such kind of MHD solutions intersect 
the magnetosonic critical curve, we are able to construct 
infinitely many such solutions and thus the collection 
of this kind of MHD solutions forms a one-dimensional 
continuum. We emphasize that this series solutions con- 
structed are closely related to the random magnetic field 
in the problem. Complementary to the hydrodynamic 
EECC solutions (Lou and Shen, 2004), the new MHD 
EECC solutions in this subsection have different asymp- 
totic behaviours at large x. The fact that solution (80 1 
also serves as an asymptotic behaviour at large x [see 
equation (83 1] should be emphasized here; we further 
confirm the validity of this asymptotic MHD solution by 
calculating one more term in the series expansion of this 
asymptotic MHD solution. 



Phase Diagram for xp = 0.2 at 7 = 1.25 and h = 1 
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Fig. 32 Density-speed phase diagram for 7 = 1.25, n = 0.75, 
h — 1 and a chosen meeting point xp = 0.2. The 'inward 
curve' is the phase path by integrating inward from different 
X {x > xf) along the magnetosonic critical curve to reach a 
chosen xp = 0.2 using the type I MHD eigensolution, and 
the 'outward curve' denotes the phase path integrating from 
different x {x < xp) along the magnetosonic critical curve 
to reach xp using the type II MHD eigensolution. There is 
no spiral pattern for the 'outward curve' and only one single 
match exists. 



4.. 3. 6 MHD Hunter Type Solutions 

Another type of MHD solutions crossing the magnetosonic 
critical curve once is here referred to as the Hunter type 
solution (Hunter, 1977) generalized to a magnetized poly- 
tropic gas and parallel to isothermal examples shown in 
figure 6 of Lou and Shen (2004). We obtain these MHD 
series solutions by solution matching in the phase dia- 
gram at a chosen meeting point x = xp where these 
two solutions have the same values of a, v, a' and v' . 
In reference to EECC solutions in the isothermal and 
unmagnetized formulation (Lou and Shen, 2004), we ex- 
pect that for 7 close to unity and h not very large, MHD 
generalizations of this Hunter type of solutions should 
exist. We explored the case of 7 = 1.01, n — 0.99 and 
h = 1, and chose a meeting point at xp — 0.5. The 



matching in the phase diagram is displayed in Figure 21 



with a familiar spiral pattern, indicating that there may 
exist infinitely many matches, corresponding to infinitely 
many discrete semi-complete MHD solutions of this type. 



4.4 Similarity Solutions Crossing the Magnetosonic 
Critical Curve Twice 

There exist semi-complete MHD solutions crossing the 
magnetosonic critical curve twice. We explore them using 
the procedure of Lou and Shen (2004) . The advantage of 
being able to reduce from three coupled nonlinear MHD 
ODEs to two will become apparent as we manage to 
match solutions in a two-dimensional phase space of v(x) 
and a{x). 

We chose different meeting points, denoted by x^, in 
the cases of 7 = 1.01 and 7 = 1.25, for a usual polytropic 
gas with n = 2 — 7. In general, higher values of h do not 
lead to matches in the range 0.01 — 0.8 for xp values 
that we have systematically searched. For 7 = 1.01 and 
smaller /i or /i = 0, we find that there are likely to be 
infinitely many matches as in Lou and Shen (2004) by 
empirical inferences, while for 7 = 1.25, we could only 
find a single match for small /i or /i = by a careful 
numerical search and by observing the variation trend of 
phase paths. 



Figure 22 presents the first three such MHD solutions Solutions with 7 — 1.01 and Small h Values 



and Figure 23 is an enlarged version of Figure 22 reveal- 
ing that the number of stagnation points increases for 
these solutions. These features show self-similar magne- 
tosonic oscillations in a magnetized polytropic gas. These 
similarity solutions are the MHD counterparts of the 
isothermal hydrodynamic solutions (Hunter, 1977; Lou 
and Shen, 2004). 



Figs. [24]through[26|present the case of 7 = 1.01, n — 0.99 
and h = Q for a nearly isothermal gas without a random 
magnetic field. At the meeting point xp — 0.5, we find 
that the 'outward curve' obtained by integrating from 
small X along the sonic critical curve to reach a chosen 



meeting point xp has a spiral pattern (Fig. 24 1 indicating 
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the trend for discrete yet infinitely many semi-complete 
solutions that cross the sonic critical curve twice. We 
also find that the number of stagnation points where a 
solution intersects the x-axis with v — increases for 
successive matches along the 'outward curve', i.e., solu- 
tion 1 has one stagnation point, and solution 2 has two 
stagnation points, and so forth (Fig. 26 1. This feature is 
very similar to the type 1— type 2 match solution in Lou 
and Shen (2004). Asymptotic behaviours of the four so- 
lutions obtained here are similar: in the small x regime, 
they approach the free-fall solutions (l50| and (51 1, while 



in the large x regime, they approach asymptotic solu- 
tions (39) and (40 1. The corresponding parameters Aq, 
Bq and m(0) are summarized in Table [s] The number 
of solutions that one can construct depends upon the 
numerical accuracy. With a numerical precision of a rel- 
ative error lO"^**, we have found at least five matches in 
a systematic and careful numerical exploration. 

Figures |27] through [29] show relevant results for an ex- 
ample of 7= 1.01, n = 0.99 and h — 1, a, nearly isother- 
mal and weakly magnetized polytropic gas, which should 
be similar to the corresponding unmagnetized isothermal 
case (Lou and Shen, 2004). Note that the third curve in 
Fig. [28] is multiplied by a factor of 5 for the compact- 
ness and clarity of the presentation. In constructing the 
phase diagram, we again find a familiar spiral pattern as 
in the h — case, and properties of stagnation number 
and asymptotic behaviours are all qualitatively similar. 
The corresponding parameters Aq, Bq and m(0) are sum- 
marized in Table [9] Figures 30 and 31 present the case 
of 7 = 1.25, n — 0.75 and h — for an unmagnetized 
polytropic gas. Within our numerical accuracy, we did 
not find a spiral phase pattern and there exists only one 
single match in this case. This match is very much like 
solution 1 in the nearly isothermal case without a ran- 
dom magnetic field. 

4-4-^ Solutions with 7 = 1.25 and Small h Values 



Figures 32 and [33| present the case of 7 = 1.25, n — 
0.75 and h — I for a weakly magnetized polytropic gas, 
which is fairly similar to the unmagnetized case oih = 0. 
The asymptotic behaviours in both cases are similar to 
previous nearly isothermal cases and the corresponding 
parameters Aq, Bq and m(0) are contained in Table 10 



4-4-3 Solutions with Larger h Values 

For a stronger magnetic field with a larger h, no solution 
match can be found in the density-speed phase diagram. 



Figures 34 and 35 are two typical phase diagrams show- 
ing the two cases of 7 = 1.01, n — 0.99, h — IQ and of 
7 = 1.25, n = 0.75 and h — 10, respectively. A possible 
interpretation of this feature is perhaps that for a large 
h, one can construct semi-complete solutions crossing the 
magnetosonic critical curve once as discussed in Section 
|4.3| while any outward solution curve will not reach the 



Table 8 The three cor responding parameters Aq, Bq and 
m(0) used in Figure 25 for hydrodynamic polytropic EECC 
solutions which are nearly isothermal. 





Values of 


7 = 1.01 and /i = 0. 


No. 


Ao 


Bo 


m(0) 


1 


5.200 


1.883 


0.3661 


2 


1.200 


-0.7958 


4.7489 X 10"'* 


3 


2.415 


0.3145 


1.096 X lO"'^ 


4 


1.904 


-0.09298 


6.6437 X 10"* 



Table 9 T he three corresponding parameters Aq, Bq and 
m(0) in Fig. 28 for MHD polytropic solutions which are nearly 
isothermal. 



Values of 7 = 1.01 and h = 1. 



No. 


Aq 


Bq 


m(0) 


1 


2.420 


0.3719 


2.441 X 10"^ 


2 


2.0155 


-0.03997 


1.304 X lO"'' 


3 


2.058 


3.868 X 10"^ 


8.590 X 10"^^ 



Table 10 T he t hree corresponding parameters Aq, Bq and 



m(0) in Figs. 31 and 33 



Specihcd parameters are 7 = 1.25 and n — 0.75. 



h 


Aq 


Bo 


m(0) 





4.399 


3.132 


2.243 X 10"^ 


1 


5.894 


2.267 


4.885 X lO""* 



magnetosonic critical curve but continue to its lower left 
to infinity to match the asymptotic MHD solution (83 1. 



5 Summary and Discussion 



We developed an MHD formulation to describe a quasi- 
spherical symmetric gas cloud obeying a conventional 
polytropic equation of state in the presence of a random 
magnetic field and searched for various possible semi- 
complete MHD similarity solutions within this model 
framework. Here, we mainly focus on MHD effects of 
a random magnetic field; other effects such as diffusions 
and radiative diagnostics will be pursued in separate pa- 
pers. 

This formulation reveals several key MHD features 
of a gas cloud when a polytropic equation of state and a 
random magnetic field are combined together. First, the 
magnetosonic critical curve diverges when x approaches 
zero, and may head up to the first quadrant in the —v 
versus x presentation when h becomes sufficiently large. 
In the absence of a random magnetic field, the sonic criti- 
cal curve also diverges as x approaches zero. In the pres- 
ence of a random magnetic field, the isothermal mag- 
netosonic critical curve remains finite at a: = (Yu 
and Lou, 2005; Yu et al., 2006). Therefore, this small 
X diverging behaviour is fundamentally related to the 
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1.5 



MHD Solutions Crossing tlie Critical Line Twice 



0.5 



-0.5 




1 



Fig. 33 The single match of the similarity MHD solution 
for 7 = 1.25, n = 0.75, h = 1 and a chosen meeting point 
Xf = 0.2, obtained with the two magnetosonic crossing points 
located at a; = 3.8496 {a = 0.2087287, v = 0.9027) and 
X = 5.6384 X 10"^ (a = 36912, v = -4.1624) along the 
magnetosonic c ritic al curve. This match is shown in the phase 
diagram of Fig. |32| 



polytropic approximation in contrast to the isothermal 
approximation. Secondly, for a sensible range of mag- 
netic field strengths, the so-called mEWCSs may de- 
scribe the large-scale MHD evolution of a magnetized gas 
cloud. But for too strong a random magnetic field, such 
a mEWCS does not exist. In this regard, our formalism 
is naturally applicable to a weakly magnetized gas cloud 
in large-scale outward and inward motions. For Bq — 0, 
a magnetized gas cloud may have a tendency to expand 
as h becomes sufficiently large. Thirdly, there exist addi- 
tional MHD EECC solutions that cross the magnetosonic 
critical curve once with asymptotic MHD solution ( 83 1 
at large x; these possible asymptotic MHD expansion 



solutions at large x were not realized before (Lou and 
Shen, 2004; Hunter, 1977). We note that semi-complete 
or shock solutions matching with such asymptotic poly- 
tropic expansion solutions may be constructed without a 
random magnetic field with h = 0. Fourthly, there exists 
a novel series asymptotic MHD solution ( 62 1 and ( 63 1 for 



small X with a very dense core as well as a strong cen- 
tral magnetic field, and one can construct semi-complete 
self-similar MHD solutions using this asymptotic MHD 
solution. Finally, for 7 = 1.25 as an example of illustra- 
tion, we have not observed the trend of spiral pattern in 
the density-speed phase diagram for solution matches. 
Therefore, there seem to be no trend of infinitely many 
semi-complete solutions crossing the critical curve twice 
under certain situations. 

Our model analysis for a magnetized polytropic gas 
cloud requires the following assumptions: the magnetic 
field is locally random on small scales with a large-scale 
quasi-spherical symmetry Zel'dovich and Novikov (1971), 



the profiles of mass density, radial velocity, gas tem- 
perature and thermal pressure are all taken to possess 
a quasi-spherical symmetry on large scales. We impose 
these assumptions as the first approximation for a quasi- 
spherical magnetized gas medium, and we try to assess 
the role of a random magnetic field under such assump- 
tions quantitatively and qualitatively. Our model analy- 
sis is specific and in details, while for astrophysical appli- 
cations, we gain a physical sense qualitatively and should 
be careful with the adopted approximation and thus the 
model limitations. An important conceptual issue should 
be noted here. Even for a hydrodynamic flow of spheri- 
cal symmetry in the absence of a random magnetic field, 
something must happen around the centre to destroy the 
self-similarity and the spherical symmetry except that 
a central black hole may accrete gas materials and ab- 
sorb them in a smooth and quiet manner (e.g., Cai and 
Shu, 2005). For any other central objects such as neu- 
tron stars, white dwarfs, or main-sequence stars, spheri- 
cal symmetric accretions will unavoidably lead to central 
activities. Therefore, a self-similar evolution of hydrody- 
namic gas flows is possible outside the influence sphere 
of central activities. Parallel to this physical rationale, 
an MHD inflow of quasi-spherical symmetry in the pres- 
ence of a random magnetic field will eventually give rise 
to a central sphere of MHD activities. Again, we expect 
that a self-similar evolution of MHD gas flows is possible 
outside this influence sphere of central MHD activities. 

In this analysis, we have adopted a polytropic equa- 
tion of state ^ . We note the existence of two different 
deflnitions for the term 'polytropic': one corresponds to 
an overall equation of state p = np'^ with a constant k 
throughout the dynamic evolution in spatial and tem- 
poral domains (Goldreich and Weber, 1980), while the 
other corresponds to the entropy conservation along the 
streamline 



d_ 
di 



d 
or 



log 



P 







(108) 



indicating that every infinitesimal portion of gas obeys 
the polytropic state equation along each streamline, but 
not necessarily in the global sense (e.g.. Bouquet et al., 
1985), i.e., the entropy per unit mass assigned to each 
streamline can be different. In the present flow problem, 
all streamlines are radial and a constant n for all stream- 
lines would certainly meet the requirement of a quasi- 
spherical symmetry. In other words, we have taken the 
former, yet we recognize that if the latter is adopted, the 
parameter n would be a free parameter, which may be 
adjusted to fit the requirements of, for examples, initial 
density profiles, initial 'equations of state' or asymptotic 
finite radial speeds at large x. By taking the former def- 
inition of a 'polytropic' gas, our analysis is primarily for 
the case of n = 2 — 7. In fact, Suto and Silk (1988) set 
the equation of state to be 



(109) 
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where K{t) [corresponding to the notation K{t) of Suto 
and Silk (1988)] takes the form of a power law in time 
t to accommodate certain unknown energetic processes. 
This form of K(t) depending on t constrains the range 
of index 7. If we adopt the usual poly tropic equation of 
state from the perspective of Suto and Silk (1988), pa- 
rameter n should take on the value of 2 — 7. If instead, k 
is allowed to assume the form of power laws in both t and 
r and equation ( 109 ) is adopted, one more free parameter 



(equivalent to our parameter n here) would emerge, as 
has been done for example by Cheng (1978) and Fatuzzo 
et al. (2004). For a careful study on behaviours of 'poly- 
tropic' gas under the latter definition, we should allow 
for the freedom of parameters in this problem and treat 
K as a function of both r and t. In short, we have treated 
here only the case of 7 = 2 — n rather than n = 1 (Suto 
and Silk, 



Goldreich and Weber (1980) studied the specific case 
of 7 = 4/3, invoking this model to describe a homologous 
evolution of the core collapse of a supernova progenitor. 
Yahil (1983) noted that even for a slight departure of 
7 value from the exact value 4/3, the core collapse will 
be no longer completely homologous. At this stage, the 
results of Goldreich and Weber (1980) are beyond the 
framework of our study, because for 7 = 4/3, our similar- 
ity MHD transformation (|9| becomes invalid. With this 
qualification, we briefly discuss limiting features of our 
model analysis as 7 approaches 4/3. In our model frame- 
work, especially in reference to nonlinear MHD ODEs 
(23 ) and (24), parameter n will approach 2/3 as 7 4/3; 
it then follows that the factor {nx — v)/{3n — 2) goes 
to infinity, except when v approaches 2x/3 simultane- 
ously, corresponding to a homologous collapse. Since first 
derivatives v' and a' in terms of x should remain finite 



in both MHD ODEs (|23| and i^, we conclude that as 
n approaches 2/3, we must require v = 2x/3 in order to 
avoid infinite first derivatives of both v and a. There- 
fore in our MHD homologous collapse should still 
exist (see Yahil 1983 for an analysis of departure from 
7 = 4/3). This argument serves only as a quick analysis 
and detailed calculations will be pursued in a separate 
paper. 

The consequences of adopting a polytropic equation 
of state and incorporating a random magnetic field in our 
formulation of gas collapses and outfiows lead to several 
novel features of possible MHD solutions in comparison 
to earlier analyses. Firstly, the new asymptotic MHD so- 
lutions (62) and (63) describe a gas cloud in which a 



random magnetic field and gas materials are both com- 
pressed in the core region with a very high total pressure 
in the core. This asymptotic MHD solution may repre- 
sent a self-similar dynamic evolution course that forms 
a compact stellar object that contains extremely intense 
magnetic field and high mass density. The formation of 
this type of objects requires an embedded random mag- 
netic field to be strong enough, because for h = this 
asymptotic solution for small x does not exist (see Yu 



et al. 2006 for this type of asymptotic MHD solutions in 
an isothermal gas medium). Secondly, with or without a 
random magnetic field, the EECC solutions crossing the 
critical curve once provides a possibility that a gas cloud 
with a collapsing core has a much faster expanding en- 
velope at a high radial speed. A more general treatment 
of fiows at relativitistic speeds under a strong gravity 
would involve the general theory of relativity; in fact, the 
expanding envelope represents a kind of nonrelativistic 
Einstein-de Sitter expansion (Shu et al., 2002; Lou and 
Shen, 2004; Yu and Lou, 2005). Here, we further spec- 
ulate the possibility that a random magnetic field and 
a polytropic equation of state might be relevant to the 
large-scale expansion of the entire universe involving a 
weak and random magnetic field. Finally, although MHD 
EECC solutions crossing the magnetosonic critical curve 
once can be constructed with a strong magnetic field, 
it can happen that the kind of EECC solutions crossing 
the critical curve twice (Lou and Shen, 2004) does not 
exist. The point is that the two MHD solutions having 
different asymptotic forms may not connect each other. 
In other words, it seems to be the case that for a very 
strong magnetic field, there are no semi-complete solu- 
tions to approach asymptotic MHD solutions ( 39 1 and 
(40) at large x. This is to say that for these cases, the 



MHD gas fiow system does not have a non-increasing 
self-similar radial fiow profile. 

We now describe and comment on the self-similar os- 
cillatory behaviour about the abscissa a;— axis described 
, of the numerical solutions Hi, 2Ii 



in subsection 4.3 
and 3Ii in Figs. 17] and [18] for the hydrodynamic and 
MHD cases, respectively. Typically, the —v versus x solu- 
tion curve behaves in the following manner. In the range 
of small x values, ~v appears to cross the x-axis repeat- 
edly and regularly for a number of infinite times, each 
cycle with a period decreasing with smaller x, and the 
vibrating amplitude also decreases as x decreases in an 
apparent power-law fashion (see Fig. 36). We explored 
this feature down to cc > 10~^° and did not find the 
curve to crash on to the singular surface. An illustrat- 
ing example of this vibration behaviour is clearly shown 
in Fig. 36 This behaviour, after careful examination and 



comparison of numerical and analytical results, turns out 
to be the type II 'quasi-static' asymptotic solution de- 
rived by Lou and Wang (2006) where the hydrodynamic 
problem in the absence of magnetic field was solved. Ap- 
pendix ]G] gives a brief account of the basic results. This 
behaviour is a characteristic of the polytropic equation of 
state with or without magnetic field. However as shown 
in Appendix ]G] the presence of a magnetic field mod- 
ifies the asymptotic solution behaviour in a non-trivial 
manner. 
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In the approximation of quasi-spherical symmetry, we take 
< ve >= and < «0 >= 0. The 9— and components of 
the magnetic induction equation can thus be simplified to the 
following forms of 



dBe 
dt 



^- — {rVrBg 

r or 



(117) 



A Magnetic Induction and the Lorentz Force 

In this appendix, we rearrange the magnetic induction equa- 
tion and discuss approximations made regarding the mag- 
netic Lorentz force in our MHD formulation with the sim- 
plification of a quasi-spherical symmetry on large scales (Yu 
and Lou, 2005; Yu et al., 2006). 

For non-relativistic quasi-neutral MHD flows, the mag- 
netic Lorentz force density can be split into two terms 



(V X B) X B 
47r 



+ 



{B ■ V)B 

47r 



(110) 



by a vector identity and the V-B = condition. The first and 
second terms on the RHS represent the magnetic pressure and 
tension force densities. The radial component of the magnetic 
tension force on the RHS is 



1 

47r 



„ dBr , Bo dBr , 

or r 06 rsm 



B^ dBr {Bl + Bl) 



(111) 



where the first term in the square brackets cancels a rel evant 
term of the magnetic pressure force in expression ( |110[ ). For 
a large-scale average of a random magnetic field, we assume 
< Bg-dBr/dr >?s^ and < B^-dBr/dO >^ 0. The nonvanish- 
ing radial component of the magnetic Lorentz force density 
associated with the mean square of the random transverse 
magnetic field is 



d_ 
dr 



< Bf > 
Stt 



47rr 



(112) 



as shown in equation 

In the absence of resistivity and ambipolar diffusions etc., 
the magnetic induction equation appears as 



dB^ 
dt 



1 d 



(rvrB^) 



r dr 

which then immediately lead to 



d ( r'^Bl\ Vr d , iryi-, , r. 



dt 



2 o2t'r 



, 



(118) 



(119) 



where Bf = Bg + B^ is proportional to the energy density 
of the random transverse magnetic field. Equation (119 1 is 
simply the magnetic induction equation ([5|. For the radial 
component of the magnetic induction equanon (Yu and Lou, 
2005; Yu et al., 2006), we have approximately 



d{r'^Br) d{r^Br) 



dt 



+ v,.- 



dr 



= r^{BfV)vr 



(120) 



Equations ( 119 1 and ( 120 1 here differ from equations (1') and 
(1) in the formulation ot Chiueh and Chou (1994). 



B Partial Derivatives of A, V and X 

The first-order partial derivatives of A, V a nd X as functions 
of a, V, and x ar e us e d in subsection |3.6[ They are derived 
from expressions (25 1, (26 1, and ( |27[ ) for A, V and X, respec- 
tively, with straightforward manipulations. We present below 
the explicit expressions of these partial derivatives, which are 
used in numerical integrations to determine the MHD solu- 
tion behaviours in the vicinity of the sonic or magnetosonic 
critical curves. 



Aa = 2a 



(n — 1)11 - 



2{x — v){nx — v) 



+3a 



2hx + 



(nx — v) 



(371 - 2) 



(121) 



dB 
~dt 



V X {~v X B) , 



(113) 



which can be written explicitly in three component forms in 
spherical polar coordinates {r,9,(f>), viz.. 



dBr _ 1 d 
dt r sin 9 dt 
1 d 



Sm9{VrBg — VgBr) 
:{v^Br — VrB^) , 



dBg 

dt 



r sm 9 
_ld_ 
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(3n - 2) 



2h 



3n- 2 



4v 
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2n + 
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(114) 14= {n- l)[v{nx - v) + Ahax'^] 

2{nx — v)'^ 2 7-1 {x — v) 



(3n - 2) 



a — 27 a 



, ,, > 2a^(nx — v) 2-ioi^ 

K = (n-l)(n.-2.)a-^A_i + JI^ 



(122) 
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Phase Diagram for xp = 0.5 at 7 = 1.01 and h = 10 
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Fig. 34 Density-speed phase diagram for 7 = 1.01, n = 0.99, 
h = 10 and a chosen meeting point xp — 0.5. The 'inward 
curve' is a phase curve by integrating inward from different 
X along the magnetosonic critical curve to xp = 0.5 and the 
'outward curve' is the phase curve integrating outward from 
smaller x to reach the same xp. There is no match in this 
phase diagram, because for small x along the magnetosonic 
critical curve, the solution curve encounters the magnetosonic 
singular surface before reaching xp — 0.5. 



Va; = (n — 1) (nv + 4/iQa;) a + 



2no?{nx — v) 2'ycPv 
(3n - 2) 
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Phase Diagram for xp = 0.2 at 7 = 1.25 and h = 10 
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Fig. 35 Density-speed phase diagram for 7 = 1.25, n = 0.75, 
h = IQ and a chosen meeting point xp — 0.2. The 'inward 
curve' is a phase curve obtained by integrating inward from 
different x along the magnetosonic critical curve to reach xp 
and the 'outward curve' is a phase curve integrating outward 
from a smaller x. There is no solution match here. 



We determine the second-order derivatives of a{x) and 
v{x) with respect to x by the following procedure. As A = 
V = X = along the sonic or magnetosonic critical curve, 
one readily obtains 



X'A" - X"A' A" ~ a'X" 



2{X') 



2X' 



(130) 



Xa = {nx — v) —70^ — 2hax 



X^ = ~2a{nx — v) , 



Xx = 2na{nx — v) ~ 2ha x . 



(127) 
(128) 
(129) 



X'V" - X"V' 



C Second-Order Derivatives of a.{x) and v{x) 

For nodal points referred to in subsection |3.6| there are in- 
finitely many solutions crossing the sonic and magnetosonic 
critical curves . Otr ly the analy tical eigensolutions determined 
by equations ( 100 1, ( 101 1 and ( 102 1 are of main interest here, 
and a numerical integration is unstable in the direction away 
from the critical curve if one uses only the first-order deriva- 
tives of a{x) and v{x) with respect to x. So one should work 
out the second-order derivatives of a{x) and v{x) with respect 
to X in order to pick out the analytic eigensolutions among 
many other solutions having weak discontinuities across the 
sonic or magnetosonic critical curves (Lazarus, 1981; Whit- 
worth and Summers, 1985; Hunter, 1986). These weak discon- 
tinuities may be regarded as weak shocks (Boily and Lynden- 
Bell, 1995). It is also possible to construct various forms of 
shocks across the sonic or magnetosonic critical curve (Tsai 
and Hsu, 1995; Shu et al., 2002; Shen and Lou, 2004; Bian and 
Lou, 2005; Yu and Lou, 2006; Lou and Wang, 2006, 2007). 



V" - v'X" 
2X' 
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We denote the second derivatives of A, V and X as 



A" 
V" 
X" 
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(131) 



(132) 



where each matrix element fij — fij{a,v,x,a' ,v') of matrix 
F is a function of a. v , x, a' and v' (i. j = 1. 2. 3). Sub- 
stituting equation (132 1 into equations (1301 and (131 1, one 
immediately obtains 

{2X' + a'f-n - /ii)a" + (q732 - fi2)v" = /13 - a'^a , (133) 



(^^731 - /2i)a" + {2X' + v'f32 - f22)v" = Hz ^ hi , (134) 

from which a" and v" can be solved directly for the non- 
degenerate case and be expressed and computed by corre- 
sponding X, a, V a' and v' at points along the magnetosonic 
critical curve. For the convenience of reference and checking, 
we summarize the explicit expressions for functions X' and 
fij below. 
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„ Vibrating Behaviour of MHD Solutions at Small x 
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Fig. 36 An example for the magnetosonic oscillation be- 
haviour of an MHD similarity solution at small x with both 
(the absolute value of v) and x displayed in logarithmic 
scales. Relevant parameters are 7 = 1.25, n = 0.75 and h — 2. 
This similarity MHD solution is obtained by integrating from 
the point a = 0.18, x = 3.7395 and v — 0.3863 on the magne- 
tosonic critical curve inward to a; = 0. The alternating plus 
and minus signs marked along the curve denote the actual 
sign of —V to compensate the loss of sign information as we 
show \v\ in a logarithmic scale here. The plot clearly suggests 
that the velocity amplitude scales with 3; in a power law. 



appropriate, we infer that SK = K{a\ + 5a, v\ + 5v, x\ -|- 
5x) — K{a\, vi, xi) = 0. It then follows from 
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with K denoting A, V or X in turn, that 
= 
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along the magnetosonic critical curve for the given set of 
(ai, v\, xi). For a nontrivial set of {5a, 5v, 5x), we must 
then require 



Ay A^ 

Vc K, 

Xa Xv X^ 



(147) 



along the magnetosonic critical cur ve, w hich is clearly equiv- 
alent to Di = in cubic equation ( 106 1. 



We conclude that because the magnetosonic critical curve is 
a continuous curve as we have solved it an alytically, it must 
lead to = in cubic equation ( 106 1 . A brief proof is 
presented belo w. Since functionals J\, V, X as defined by 
equations l\25\ — ( |27[ ) together determine the magnetosonic 
critical curve (i.e., setting A = V = X — simultaneously), 
we begin by assuming that the point (ai, vi, xi) in the 
a — v — x space is on the magnetosonic critical curve such that 
7^(01, vi, xi) = where K denotes A, V or X, respectively. 
We further assume that the set (qi + 5a ,vi + 5v , a;i -|- 
5x) is also along the magnetosonic critical curve and is very 
near to the former point. Then from the requirement K{ai + 
5a ,vi + 5v ,Xi + 5x) = 0, with K denoting A, V and X as 



E Determination of Magnetosonic Critical 
Curves 

One can numerically determine ranges of a values for which 
the signs of yli, Bi and Ci and of the determinant A = Bf — 
4:Ai Ci in equation ( 35 1 do not change and one can thus decide 



which root of x^ should be picked up in equation ( 34 1 and how 
many roots are physically relevant once specific values of 7, h 
and n are prescribed. This numerical exploration is important 
as one needs to know the sonic or magnetosonic critical curve 
in the entire range of 0"'" < q < +00. We have analyzed the 
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Table 11 Algebraic signs for quantities Ai, Bi, Ci and A, 
and the corresponding root numbers for a usual polytropic 
gas with 7 = 1.01 and n = 0.99 



Table 12 Algebraic signs for quantities Ai, Bi, Ci and A, 
and the corresponding root numbers for a polytropic gas with 
7 = 1.25 and n = 0.75 



7 


= 1.01, n = 0.99, /i = 


7 


= 1.25 , n = 0.75 , h = 


ail 


= 0.202728 , ai2 = 0.282682 


asi 


= 0.279006 , 052 = 0.368686 


a G 


(0 , an) (an , ai2) (ai2 , +co) 


a G 


(0 , asi) (asi , 052) (ass , +00) 


signs 


Bi>0 Bi, A>0 


signs 


i^i > iii , Z\ > 




all < others < Ai, Ci < 




all < others < Ai , Ci < 


root No. 


2 


root No. 


2 



7=1.01, n = 0.99, h = OA 
aai = 0.048500, 022 = 0.103907, 023 = 0.202728 




7 = 1.25, n = 0.75, h = l 
= 0.09283 , a62 = 0.174661 , aes = 


0.279006 


aG (0, a2i) (a2i, 022) (022, 023) (023, +00) 


a G 


(0, aei) (asi, 0^2) (a62, aes) 


(a63, +00) 


signs A>0 Ai, A>0 Ai, Bi, A > 
all < others < Bi, Ci < Ci < 
root No. 1 1 


signs 
root No. 


A>0 Ai,A>0 
aU < others < Bi, Ci < 
1 


Ai,Bi,A > 
Ci < 
1 





7 = 1.01, 
aai = 4.9485 x 10" 



n = 0.99, 
a32 = 0.05 1866, 



a33 = 0.202728 



a G 


(0, aai) (asi, 032) (032, 033) (033, +<x) 


signs 
root No. 


A>0 Ai,A>OAi,Bi,A>0 
all < others < Bi, Ci < Ci < 
1 1 




7=1.01, n = 0.99, h = W 
a4i = 4.9440 x lO""*, 042 = 0.027885, 043 = 0.202728 


a G 


(0, a4i) (041, 042} (a42, 043) (a43, +00) 


signs 
root No. 


A>0 Ai,A>0 Ai,Bi,A>0 
all < others < Bi, Ci < Ci < 
1 1 



usual polytropic cases of 7 = 1.01, n = 0.99 and 7 = 1.25, 
n = 0.75 for different h values. Ifere, we summarize the main 
result s. T he case that 7 = 1.01 and n = 0.99 is contained in 
Table [TT]and the case that 7 = 1.25 and n = 0.75 is contained 
in Tablell2l 



F MHD Eigensolution Behaviours in the 
Vicinity of the Magnetosonic Critical Curve 

Similarity MHD solution behaviours in the vicinity of the 
magnetosonic critical curve can be determined by the signs 
of both A = B ^ — 4 C3 and C3 as defined by equation ( 106 1 
(see subsection 



3.6 for details). One can numerically deter 
mine the values ot A and C3 given specific x values of dif- 
ferent magnetosonic critical curves corresponding to different 
parameter sets of 7, n and h. We describe the main results 
of this exploration for the cases of 7 = 1.01, n = 0.99 and 
7 = 1.25, n = 0.75 with different h values. In Table [131 
present the relevant results for 7 = 1.01 and n = 0.99, while 
in Table [141 we sho w th e results for 7 = 1.25 and n = 0.75. In 
both Tables [13] and [lH the sign 's' indicates a saddle, 'n' indi- 
cates a node, and 'c' indicates a centre or a spiral regarded as 
being unphysical (Jordan and Smith, 1977). Note in particu- 
lar that among the cases we have explored and not displayed 
in the above Tables, there exist such magnetosonic critical 
curves that all points along the lines are saddle points. These 
include the cases of 7 = 1.01, n = 0.99, h — 10 and 100, and 
of 7 = 1.25, n = 0.75, h — 1, 10, and 100. There are also cases 
of /i > with the coexistence of nodal, saddle and spiral or 
centre points along the magnetosonic curve, e.g., the case of 
7 = 1.01, n = 0.99 and h = 0.0001. The fact that there exist 



7 = 1.25, n = 0.75 
an = 0.0090897 , 072 = 0.107674 

"a (0 ,071) 



signs 
root No. 



all < 




(a7i , 072) 

zi > 
others < 




a73 = 

(a72 , aTsY 
Ai,A>0 
Bi , Ci < 



0.279006 

(073 , +00) 



Ai,Bi,A> 
Ci < 
1 





7 = 1.25, n = 0.75. 


h = 100 


asi = 8.9955 x lO"* , 


ias2 = 7.1799 x 10"^ , ass = 0.038660 


a G 


(0 , agi) 


[asi , Q82) 


(a82 , aga) 


signs 




A>0 


Ai , A>0 




all < 


others< 


Bi , Ci < 


root No. 








1 




084 = 0.079160 , as5 = 


0.279006 


a G 


(ags , 034) 


(a84 , O85) 


(085 , +00) 


signs 


A>0 


Ai , A>0 


Ai, Bi, A>0 




others < 


Bi , Ci < 


Ci < 


root No. 





1 


1 



spiral or centre points in unmagnetized cases should be noted 
because there would be no global semi -com plete smooth so- 
lution crossing such points (subsection 3.6 1. 



G Asymptotic Behaviours Approaching the 
Quasi Magnetostatic Solution 



As mentioned above, global magnetostatic solution (741 also 
serves as an asymptotic behaviour for small x. This asymp- 
totic solution turns out to be the only solution characterized 
by a n approaching zero faster than 0{x) and o ~ x~'^^" 
when n = 2 — 7 for a usual polytropic gas. In order to deter- 
mine the behaviour of this asymptotic solution, we assume to 
the leading order 



V = Lx^ + ■ 



(148) 



(2-7) 



1 1/(7-2) 



^-2/(2-^)^^^^, 



27(4-37) 7 

where the next-order a variation Aa is given by 
Aa = iVx^-^-^/" , 



(149) 



and L, K and A'^ are three constant complex coefficients. 
Apparently, Re{K) > 1 is required for v to approach zero 
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Table 13 Algebraic s igns of both determinant A = Bf — 4C3 
and C3 for equation ( 106 1 and type identifications for local 
behaviours of MHD eigensolutions in the vicinity of magne- 
tosonic critical curves for a conventional polytropic gas with 
7 = 1.01 and n = 0.99, where 's', 'n' and 'c' denote a sad- 
dle, node, and a centre (or spiral), respectively (Jordan and 
Smith, 1977). 



7 = 1 


01, n = 


0.99, 


/i = 




Xii = 


-- 1.0455 , 


X12 = 


= 47.31 




X G 


(0 , xii) 


(a;ii 


, 2:12) 


(2:12 , +co) 


signs 


A>0 


^ , 


C3 > 


C3 >o 




C3 <0 






A<0 


type identification 


s 




n 


c 



7 = 1.01, n 


= 0.99, h 


= 0.1 


X21 = 


1.1497 




X G 


(0 , 121) 


(a;2i , +00) 


signs 


A>0 


A, a>o 




C3 < 




type identification 


s 


n 




7 = 1.01, n 


= 0.99, 


h = 1 


X31 = 


= 1.9182 




X e 


(0 , X31) 


(xai , +00) 


signs 


A>0 


A, C3>0 




C3 <0 




type identification 


s 


n 



Table 14 Algebraic signs of b oth determinant A = Bf — 
4C3 and C3 for equation 1 106 1 and type identifications for 



local behaviours of MHD eigensolutions for a conventional 
polytropic gas in the vicinity of magnetosonic critical curves 
as 7 = 1.25 and n = 0.75, where 's', 'n', and 'c' denote saddle, 
node, and centre (or spiral), respectively (Jordan and Smith, 
1977). 



among the three complex constants K, L and A''. Clearly, K 
satisfies the following quadratic equation 

f{K) = \n 12 + n(3n - 2)ft]A"^ 
-(4 - 3n)[n/2 + (3n - 2)h\K 



+n + 7(2/n - 2)(3n - 2)/i = 



(151) 



Once a proper root of K is chosen, the complex ratio of L to 
N is determined accordingly. We emphasize that the values 
of A'" and L are not determined a priori, i.e., L can take on 
any reasonable value, while A^ is determined by their ratio or 
vice versa. 

We introduce the handy notation 



ho 



[(3 + 2V2)n - 4] [4 - (3 - 2V2)n] 



(l-n)(-. 



2n(3n - 2) 
1^ + 24n - 16) 



he 



(152) 



and reach the following conclusions. When 12 — 8\/2 < n < 
0.8 and ho < h < he, or when 2/3 < n < 12 — 8^2 for 
arbitrary h values, there exist two real roots of f{K), both 
are larger than unity; these are referred to as type I 'quasi- 
static' asymptotic solutions. When 12 — 8y/2 < n < 0.8 and 
h < ho, there exists a complex root K with its real part 
larger than unity; this is referred to as type II 'quasi-static' 
asymptotic solution. For a complex K = 7^1-1-17^2, we simply 
have 



= Lx ^ exp{iK2 Inx) 



(153) 



In the limit of /i = 0, the above MHD results reduce to those 
of a hydrodynamic analysis (Lou and Wang, 2006). 

When n 7^ 2 — 7 for an unconventional polytropic gas, an 
asymptotic similarity MHD solution in the form of i; ~ 0{x) 
and Q ~ x~^^^^~~''' also exists, viz. 

2(2 - 7- n) 
(4 - 37) 



7 = 1.25, n = 


0.75, /i = 




X41 = 11.724 , 


X42 = 11.724 




X (0 , X41) 


(^41 , 2:42) 


(a;42 , +00) 


signs A > 


Zi , 63 > 


C'3 > 


C3 < 




A<0 


type identification s 


n 


c 



(1-7) 



h + 



(2-7) 
7(3n - 2) 



X ^-2/(2-7)+, 
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-1/(2-7) 



(154) 
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faster than ~ 0{x) and Aa being the next-order term of this 
asymptotic solution as x approaches zero. 



Substituting asymptotic expr essi ons ( 148 \ and ( 149 \ into 
nonli nea r MHD ODEs (111 and ( 16 1 and using integral rela- 



tion (21 1, we obtain the tollowing two relations 



n'(l-f A-) 
2(3n- 2) 



-|- nhK — 7/1 
(1-7) 



N 



2j{3n - 2) 



n{K -1)N 
(A'-2-2/n) 



+ 



7 



-l/r. 



(3n-2) ' 



27(3n - 2) 



-l/n 



(150) 
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